arXiv:1508.07465vl [cond-mat.soft] 29 Aug 2015 


Review: Dynamics of Crowded Macromolecules/ 
Interacting Brownian Particles 

George D. J. PhillieJ^ 

Department of Physics, Worcester Polytechnic Institute, Worcester, MA 01609 

Abstract 

I review theoretical treatments of diffusion in crowded (i. e., non-dilute) solutions of globular 
macromolecules. The focus is on the classical statistico-mechanical literature, much of which 
dates to before 1990. Classes of theoretical models include continuum treatments, correlation 
function descriptions, generalized Langevin equation descriptions, Smoluchowski and Mori-Zwanzig 
descriptions, and a brief but encouraging comparison with experimental results. The primary 
emphasis is on measurements made with quasi-elastic light scattering spectroscopy; I also discuss 
outcomes from fluorescence photobleaching recovery, fluorescence correlation spectroscopy, pulsed- 
gradient spin-echo nuclear magnetic resonance, and raster image correlation spectroscopy. I close 
with a list of theoretical papers on the general topic. 

This manuscript began as a chapter for someone else’s book. For reasons not relevant here, the 
following will not see print there. I am therefore making a distribution to interested parties. Given 
the size, this is not a journal article; I am open to publication offers. In the long run, this article 
may become a chapter in my volume Theory of Polymer Solution Dynamics, now in preliminary 
stages. Reports of typographic errors, points where the discussion is obscure, additions to Appendix 
C, and requests for additions are collegially welcomed. 
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INTRODUCTION 


A. Plan of the Work 


Recently, there has been increased interest in crowding, the effect of intermacromolec- 
ular interactions on the dynamics of macromolecules in non-dilute solutions. There does 
not appear to be a hrm contact between modern studies of diffusive motion in non-dilute 
macromolecule solutions and the classical theoretical and experimental literature on this 
topic. There have, after all, been more than four decades of intensive theoreticalll-3| and 

fin f 

experimentally, l5[ studies on diffusion by non-dilute macromolecules, and by dilute probes[6- 
9|] in nondilute solutions of proteins or long-chain random-coil polymers. 

This article focuses on one aspect of the historic literature, namely theoretical studies of 
diffusion by interacting macromolecules, especially as studied by quasielastic light scattering 
spectroscopy, fluorescence correlation spectroscopy, pulsed-gradient spin-echo NMR, raster 
image correlation spectroscopy, and related techniques. In particular, I treat the mutual, 
self, and probe diffusion coefficients Dg, and Dp of hard spheres. A few results on 
bidisperse and charged systems are noted. For interacting particles, these diffusion coeffi¬ 
cients depend on the concentration of the diffusing solutes. The concentration dependences 
reflect a complicated interplay of many effects, including inter-particle direct and hydro- 
dynamic interactions, reference frame corrections, and correlations between Brownian and 


3 
















driven motions. Some comparisons are made with experiment. 

Calculations of diffusion coefficients may be categorized by their general approach. The 
discussion in this article is partitioned by those categories. Classical treatments are macro¬ 
scopic, use the local concentration c(r, t) as a primary variable, and treat diffusive ffuxes as 
being driven by so-called “thermodynamic forces” and hindered by “dissipative coefficients”. 
The several microscopic approaches treat diffusing macromolecules as individual particles 
whose motions are modihed by their intermacromolecular interactions. Included among the 
microscopic approaches are correlation function descriptions, generalized Langevin equation 
descriptions, and calculations based on the Mori-Zwanzig and Smoluchowski equations. 

In the remainder of this section, we describe the types of diffusion coefficient. We then 
consider light scattering spectroscopy and methods for interpreting light scattering spectra. 
Section II treats continuum models for diffusion, including two- and three-component solu¬ 
tions and reference frame corrections. Section III treats correlation function descriptions, 
which are the starting point for microscopic calculations. Section IV uses Langevin-type 
equations to evaluate interesting correlation functions, including careful attention to the 
subtle correlations between the Brownian and direct components of the force on each parti¬ 
cle. Section V briefly treats Smoluchowski and Mori-Zwanzig type calculations. Section VI 
is a discussion, including a short comparison with experimental tests. 

B. Diffusion Coefficients 

The description of diffusion in terms of diffusion coefficients arises from classical ex¬ 
periments that observe diffusion over times that are very long relative to all microscopic 
molecular processes in solution. With light scattering spectroscopy, one can observe diffu¬ 
sive processes over short times and small distances, in which case information about diffusion 
more detailed than that given by the diffusion coefficient may be obtained. 

Operationally, one can identify at least three different translational diffusion coefficients. 
The mutual (or inter-) diffusion coefficient Dm characterizes the relaxation of a concentra¬ 
tion gradient. The self or tracer diffusion coefficient Dg describes the motion of a single 
macromolecule through a solution containing other macromolecules of the same species. 
The probe diffusion coefficient Dp determines the diffusion of an identihable, dilute species 
through a complex fluid. There is also a rotational diffusion coefficient that character- 
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izes whole-body reorientation; this diffusion coefficient is accessible via depolarized light 
scattering}^. Diffusion in solutions containing more than one solute component requires 
cross-diffusion coefficients for a complete characterization. 

Quasi-elastic light scattering spectroscopy (QELSS) of monodisperse solutions measures 
Dm- Dm can also be measured with a classical diffusion apparatus in which the disappearance 
of a macroscopic, artihcially induced concentration gradient is observed. The needed macro¬ 
scopic, artihcial concentration gradient can be produced in an analytical ultracentrifuge, 
permitting the mutual diffusion coefficient to be measured during a sedimentation-diffusion 
experiment. A true self-diffusion coefficient cannot be measured with QELSS, because laser 
light scattering is a coherent process. 

Light scattering may be made incoherent via studying inelastic re-emission by fluorescent 
groups, either floating freely in solution or covalently bonded to larger molecules of inter¬ 
est. Fluorescence Photobleaching Recovery (FPR), Fluorescence Correlation Spectroscopy 
(FCS), and Raster Image Correlation Spectroscopy (RICS) take take advantage of this in¬ 
coherent scattering to measure Dg of fluorophores and fluorescently-tagged particles. Alter¬ 
natively, one may measure Dg by resorting to Pulsed-Gradient Spin-Echo Nuclear Magnetic 
Resonance (PGSE NMR). 

Some techniques (e. g., fluorescence correlation spectroscopy) determine the diffusion co¬ 
efficient of a labelled macromolecular species. If the system under study contains both dilute 
labelled macromolecules and also unlabeled macromolecules of the same species (perhaps at 
an elevated concentration c), and if the label does not perturb macromolecular motion, the 
apparent Dm approaches closely Dg of the macromolecules at the concentration c 12|. 

There is no physical requirement that the labelled (“probe”) and unlabelled (“matrix”) 
macromolecules must except for the label be the same. If the probe and matrix species 
are different, one says that one is studying probe diffusion. If any of several stratagems for 
separating scattering due to probe particles from scattering due to the matrix solution is 
effective, QELSS may be used to measure the diffusion of probe particles through a back¬ 
ground matrix solution. For example, if the matrix molecules match the index of refraction 
n of the solvent, so that the matrix molecules’ ^ vanishes, and if the probe is dilute, QELSS 
is readily used to measure Dp, the self-diffusion coefficient of the probes. In other cases, 
subtraction - at the level of the held correlation functions - of spectra of the matrix solu¬ 
tion from spectra of matrix:probe solutions has permitted isolation and interpretation of the 
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spectra of diffusing probes [1^. 


C. Quasi-Elastic Light Scattering Spectroscopy 


Quasi-elastic light scattering spectroscopy, including both theoretical issues and exper¬ 
imental considerations, has been the subject of a series of monographs, including volumes 

Crosiganni, et al.[^, Cummins and Pike 18, 19|. 


from Berne and Pecorall4 
Pecora 


Chu 


Q,Q, 


20|, and Schmitz j2l(|. Readers are referred to these volumes for extended treatments 


of how the technique works. I present here only a very short summary. 

Experimentally, in a QELSS system the liquid of interest is illuminated with a laser 
beam. A series of lenses and/or irises is then used to collect the light scattered by the liquid 
through a narrow range of angles. The intensity I{t) of the scattered light fluctuates. The 
intensity fluctuations are monitored using a photodetector and photon counting electronics. 
The actual signal being analysed is the count n* of photons observed in each of a series of 
time intervals {t, t + 6t). 

The time-dependent intensity I{t) is used to determine the intensity-intensity time cor¬ 
relation function 

C{T) = {I{t)I{t + T)). (1) 


Here (■ ■ •) denotes an averaging process. 

The information about the liquid appears in the time dependence of C{t). Intensity 
measurements are made by photon counting; a real digital correlator actually determines 

all 

C{t) = ^ UiTli+r (2) 

i=l 

where Uj and are the number of photons that were counted in time intervals here 
labelled i and i + t, and where the sum is over a large number of pairs of times i and i + r 
separated by the delay time r. 

Scattering from a solution of Brownian particles is said to arise as scattering from a 
series of scattering centers. For simple spherical particles, one has a scattering center at 
the center of each particle. For polymers, which are not treated in this paper, the single 
scattering center is replaced with a line of scattering centers located along each polymer 
chain. The light that was scattered in the right direction is approximated as proceeding, 
without being scattered again, to the detector. This approximation is the first-order Born 
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approximation for scattering. There is an extensive theoretical treatment, rarely invoked 
for QELSS, for scattering beyond this simple but usually adequate approximation; for a 
systematic treatment, see Kerker 22], 

As is shown in the standard sources, the fluctuating intensity I{t) and its time correlation 
function C{t) are determined by the locations of the scattering particles via the dynamic 
structure factor 

S{q,T) = ( ^ aiajakaiexp{iq-[ri{t)-rk{t)+ri{t + T)-rj{t + T)])l . (3) 

In this equation, each variable in the quadruple sum proceeds over all N particles, q is the 
scattering vector, af is a scattering cross-section including all constants needed to convert 
from particle positions to intensities, and rj(t) and rj{t + t) are locations of particles i and 
j at times t and t + t, respectively. 

The actual correlator output is a bit more complicated than is suggested by this equation. 
In a simple linear correlator, in which the rii are all collected over equal time intervals, the 
actual output C'(r) of a digital auto correlator is a wedge-weighted average of 5'(g, r) over a 
range of delay times, the center of the wedge being the nominal correlation channel location 
and the width at half-height of the wedge being the time spacing Ar between correlator 
channels. As an exception, (7(0) is not a good approximation to S'(g, 0). In a modern multi- 
tau correlator, with increasing r the rij are aggregated into counts covering longer and longer 
time intervals. In a multitau correlator the averaging over delay times is more complicated. 
The delay time r to be assigned to a channel requires careful analysis 23| . With respect 


to the importance of careful analysis, for which see ref. 


23l. suffice it to say that the so- 


called ’half-channel correction’ is totally wrong for linear correlators, and at best a crude 
approximation for multitau correlators. Identifying C'(r) with S{q,T) is possible, so long as 
the two values of r are properly adjusted. This averaging matter is a purely experimental 
issue; in the following we assume that it has been handled correctly. 

Closely related to the dynamic structure factor is the held correlation function or inter¬ 
mediate structure factor 

g^^\q, r) = exp(fq ■ [ri(f) - + r)])^ (4) 

Under normal experimental conditions (this restriction does not include some modern mi¬ 
croscopic techniques) all information in S'(g, r) is actually contained in g^^\q,T), so it is 
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with the evaluation of the much simpler in various guises that we are concerned 

in this paper. This article conhnes itself to scattering from large-volume systems, for which 
g^^\qiT) can be extracted from S{q,T) via 

S{q,T) = A{g^^\q,T)f + B. (5) 


A is an instrumental constant. B is the baseline, the numerical value of S{q,oo) to which 
the correlation function decays. 

The transition in eq. [5] from S'(g, r) to g^^\q,T) is an extremely good approximation. 
Crosignani, et al.[l7| provides the justification. Terms of the quadruple sum of eq. [3] are 
only non-zero if the particles are close enough for their positions or displacements to be 
correlated. Terms that put one or three particles in a scattering volume average to zero. 
So long as the regions over which particle positions and displacements are correlated - the 
correlation volume - are much smaller than the region from which scattered light is being 
collected, the sum has many more terms that put two particles in each of two correlation 
volumes than it has terms that put four particles in a single correlation volume. The error 
in the approximation in eq. [5] is in the mistreatment of terms putting four particles in the 
same correlation volume; this error is small. 

The r —?• 0 limit of the pair correlation function is the static structure factor 

INN \ 

g^^\q) = g^^\q.^) = 5^ aiajexp(iq- [ri(t) - rj(t)]) \ , (6) 

\ i=l / 

which describes equal-time correlations between the positions of particles in the solution. 
In eq. El the first and second summations are, respectively, the self and distinct terms of 

c/d)(g). 


D. Spectral Analysis 


How does one extract experimental parameters from S'(g,f)? The first and most impor¬ 
tant issue is that light scattering spectra are relatively featureless, and therefore the number 
of parameters that can be extracted from a spectrum is extremely limited 2^ . From one 
spectrum, a half-dozen parameters is often optimistic. A spectrum having well-separated 
relaxations spread over four orders of magnitude in time and a very high signal-to-noise 
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ratio may yield as many as eight parameters, though the scatter in repeated measurements 
will be substantial. 

Four noteworthy approaches to analyzing S{q, r) are cumulant analysis, lineshape htting, 
moment analysis, and inverse Laplace transformation. We consider these seriatim. Cumu¬ 
lant analysis begins with the observation that the held correlation function can be written 
formally as a sum of exponentials 


POO 

j«(«.t)=/ <;r/i(r)exp(-rT) (7) 

Jo 

Eq. [7] is the statement that g^^\q,T) can be represented as a Laplace transform; it has no 
physical content. The relaxation distribution x4(r) is normalized so that 



dTA{T) 


1 . 


The hrst moment of x4(r) is 


The central moments of ^(T) are 


POO 

f = / dTAiT)T. 

Jo 


( 8 ) 

(9) 




dTA{T){T - T)^. 


( 10 ) 


In cumulant analysis, as introduced by Koppel 
the spectrum is represented by 


25[ | and greatly improved by Frisken 26j , 


( 11 ) 

\n=0 / 

The Kn are the cumulants, which are not the same as the central moments except for small 
n, Ko being the spectral amplitude and Ki = T being the average decay rate. 

A traditional advantage of the cumulant expansion is that it can be applied by making a 
weighted linear least-squares ht, namely 


ln(^(i)(g, r)) = ln{{S{q, r) - 5))^/^ = ^ - (12) 

n! 

n=0 

The In and square root change the statistical weights to be assigned by the spectral htting 
program to diherent data points. is a truncation parameter, the order of the ht; it hxes 
the highest-order cumulant to be included in the htting process. 
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Frisken 


26[ | provides an alternative expansion 


S{q,T) = 


exp(-Kit) 1 + + 


-\ 2 


+ B 


( 13 ) 


2! 3! 

in which the /i„ are the central moments. In this expansion, Ki and the fin) are to be 
extracted from S{q, t) via non-linear least sqare fitting methods, e.g., the simplex algorithm. 
(There are two very different classes of simplex algorithm. The useful one here is the Nelder- 
Meadj^ functional minimization approach. The other simplex approach is Dantzig’s linear 


programming method 28|.) The first three Kn are the same as the first three fin, but the 


higher-order Kn and fin are different. 

Before the series are truncated, equations fTTIfT^ are all ways to write g^^\q,T) as a 
convergent exact power series, using expansion coefficients that are readily described in terms 
of the Laplace transform A(r) of g^^"^{q,T). The inhnite series are therefore appropriate to 
describe an arbitrary ^(T). Some truncated series may have numerical stability issues if r 
covers a wide range of times; the Frisken expansion of eq. [I3] avoids these. Under practical 
conditions, one fits the complete g^^\q,T) or S'(g,r) to eq. [TTl eq. [121 or OQ- US] while 


varying N upwards from 1, and uses the smallest N such that further increases in N do not 
substantially improve the quality of the £t. 

The mutual diffusion coefficient obtained by light scattering is defined to be 


Dn, = (14) 

Cumulants can be written as the logarithmic derivatives of g^^\q, t) in the limit of small 
time delays, so that formally 

\ n 

H9^^\qK)]- (15) 

The presence of the limit in front of the derivative sometimes leads to the false assertion that 
cumulants only capture fast relaxations. Literal application of this formula in theoretical 
calculations of the Kn is complicated by time scale issues. The physically observable g^^^ (g, r) 
is only obtained for r greater than the correlator channel width At, so the small-time limit 
in eq. da does not include phenomena that relax completely in times At. 

An alternative to cumulants analysis is a forced fit of an assumed form to S'(g, r) or 
g^^\q, t). In the earliest days of the field S{q, r) was often force-fit to a single exponential. In 
complex fluids and glassy liquids, relaxations often take the form of a stretched exponential 
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exp(—ar^) in time, or a sum of several stretched exponentials in timejw, 29|. For interacting 
systems, the long-time part of the spectrum is sometimes found to follow a power-law decay 

^(g,r) = ar^, (16) 

a and z/ being htting parameters. Long-time power-law decays plausibly arise from 
mode-coupling behavior. Power-law tails have been observed experimentally for light 
scattering spectra of strongly-interacting charged polystyrene spheres at very low salt 


concentrations. 


3 ^. 


For some purposes it is useful to determine the moments of the held correlation 
function, where for n > 1 

0 ) • 

The moments have the physical advantage that they are dehned as integrals, not polynomial 
hts or derivatives, so that (if a functional ht can be used to extrapolate r) to r —)■ oo) 

moments and cumulants are susceptible to different sorts of noise. 

Moments give average diffusion coefficients, with the slowest decays weighted most heavily 
in the average. If the spectrum is written t) = f (iFy4(F) exp(—Fr), with A(F) having 

been normalized so that f (iF74(F) = 1, then Mi = (1/F), and similarly for the higher-order 
moments. 

In principal, spectra of interacting systems can also be ht to sums of exponentials via 
inverse Laplace transform methods. However, inverse Laplace transforms are ill-posed, so 
the outcomes are highly sensitive to noise. Some software reports the smoothest (in some 
sense) function that is consistent with the observed spectrum. The relationship between 
the smoothest function and the actual Laplace Transform can be obscure. Furthermore, 
most theories do not naturally lead to multi-exponential forms for spectra of monodisperse 
interacting systems. 

We now advance to treatments of dihusion, beginning with the continuum treatments. 


II. CONTINUUM TREATMENT OF LIGHT SCATTERING 

A. Introduction to Continuum Treatments 

Continuum descriptions of concentration may be traced back to Fick’s original memoir on 
dihusion. Continuum descriptions have been used extensively to analyse dihusion problems. 
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and are sketched here. Continuum models generally omit details of the interactions between 
individual particles. 

Many continuum models are based on non-equilibrium thermodynamics. Non-equilibrium 
thermodynamics asserts that diffusion currents may be written as products of forces and 
dissipation coefficients. This level of refinement was apparently adequate for early 20*^ cen¬ 
tury treatments of sedimentation and electrophoresis, in which particles moved under the 
influence of an external field that was independent of molecular coordinates, while dissi¬ 
pation coefficients were described in terms of averages over molecular coordinates. Some 
non-equilibrium thermodynamic models postulate non-mechanical “thermodynamic driving 
forces”, whose magnitudes are not computed from classical statistical mechanics as applied 
to microscopic molecular systems. 

Non-equilibrium thermodynamics is not consistent with the microscopic treatments given 
in later sections. If both the force and the dissipation are determined by molecular coordi¬ 
nates, then the current should be given by an average over the product of an instantaneous 
force and an instantaneous mobility, not by the product of an averaged force and an aver¬ 
aged mobility coefficient. Nonetheless, continuum models are invoked frequently. Having 
provided a necessary caveat for the following discussion, we turn to continuum models. 


B. Formulation of Continuum Treatments; Einstein Model 


In continuum treatments of diffusion, a solution is treated as having a continuum con¬ 
centration c{r,t). In the continuum treatments, r refers to a location in space, while c{r,t) 
refers to the density of scattering particles near r at time t. This description has been 
extensively used to analyze diffusion problems. Note that in this Section the coordinate r 
has entirely changed its meaning. In the previous section, the N scattering particles were 
assigned time-dependent coordinates {ri,...,rAr} that specified their locations as functions 
of time. In this section, each r is a fixed location in space. In older works on diffusion and 
experimental methods for its study, authors will sometimes jump back and forth between 
these two meanings of r. 

A simple continuum treatment of diffusion uses the continuity equation 


(9c(r, t) 
dt 




( 18 ) 
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and Pick’s diffusion equation 

J(r, t) = —Zi)i?Vc(r, t) (19) 

to obtain 

^ = ( 20 ) 

Here J(r,t) is the time- and position-dependent concentration current and Dp is the Pick’s 
Law diffusion coefficient. In eq. [20l terms in 'VDp = have been omitted, which is 

appropriate if the concentration dependence of is small enough. 

The simplest continuum model for diffusion is due to Einstein. Despite its simplicity, the 
model is adequate to predict a diffusion coefficient for dilute particles in terms of mechanical 
properties of the solute and standard statistico-mechanical considerations. Einstein consid¬ 
ered a solution of highly dilute particles, each having buoyant mass m, floating in solution. 
There is a gravitational held, strength g, so the potential energy of each particle is 

U = mgz, (21) 

with z being the vertical coordinate. Prom standard statistical mechanics, the equilibrium 
concentration of particles in solution is 

c{z,t) = Coexp{—(3mgz), (22) 

in which cq is the concentration of particles at z = 0, which is presumed to be inside the 
solution. Here (3 = where ks is Boltzmann’s constant and T is the absolute 

temperature. 

However, from Pick’s Law there will be a diffusion current 


JzD{z,t) = DF/3mgcoexp{-/3mgz), 


(23) 


the current being upward (for rUg > 0, as is not always the case) because the concentration 
is larger as one goes farther below the surface of the liquid. As a result of the gravitational 
force mg, the solute particles fall at their terminal velocity mg/f, leading to an induced 
sedimentation current 


Jzs{z,t) = ^coexp{-Pmgz), 


(24) 


In this equation / is the drag coefficient for motion at constant velocity. Even though each 
particle falls at the same terminal velocity, the sedimentation current depends on because 
the concentration of particles depends on ; 2 . 
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We are at equilibrium, so the concentration does not depend on time, meaning the diffu¬ 
sion and sedimentation currents must add to zero. The only way the two currents can add 
to zero is if 

Df = (25) 

Equation [25] is the Einstein diffusion equation. When combined with Stokes’ Law 


/ = Gnrja 


(26) 


for the drag coefficient of a sphere, one obtains the Stokes-Einstein equation 

ksT 


Dj.= 


Qurja 


(27) 


for the diffusion coefficient of a sphere. Here rj is the solution viscosity and a is the sphere 
radius. Observe that the derivation is a ’just-so’ story. The value of Dp was not calculated 
directly, for example from Newtonian mechanics. Instead, it was shown that Dp had to 
have a certain value, or the laws of classical mechanics and statistical mechanics would be 
violated. 

How is the scattering intensity related to the concentration of particles in solution? The 
key step is to recognize that the scattered light is determined by the q*^ spatial Fourier 
component aq{T) of the solution concentration, namely 


N 

= exp(iq ■ ri(r)) (28) 

i=l 

so if there is only one solute component the intermediate structure factor becomes 

g^^\q,T) ~ a^(ag(0)a_g(r)). (29) 

Here is proportional to the particle scattering cross-section, includes all the constants, 
some substance-dependent, that show how the scattered light intensity depends on the size 
of the concentration fluctuations. 

The Fourier components satisfy 

c(r, t) = Y1 ^ (30) 

g 

(pq being a time-dependent phase unique to each Fourier component. The phase disappears 
if: The aq{t) are made complex, the cosine is replaced with a complex exponential, and the 
real part is taken. 
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Substitution of eq. [30] into eq. [20] gives the temporal evolution of the as 

aq{t) = aq{0) exp{—DFq‘^t). (31) 

The scattered held mirrors the behavior of the concentration huctuations via 

g^^\q,T) ~ a2([ag(0)]2)exp(-T)j.gV). (32) 

Comparison with the models below identihes Dp as the mutual diffusion coefficient Dm- 


C. Three-Component Systems 


By extending this treatment to a 


hree-component solute:solute:solvent mixture, self- and 


probe diffusion can be discussed { 2 , Denoting the solute components as A and B with 

concentrations ca and cp, respectively, the corresponding scattering cross-sections as a a 
and as, and the amplitudes of the corresponding spatial Fourier components of the two 
concentrations as aqAit) and aqsit), respectively, the electric held of the scattered light can 
be written 

Es{q, t) ~ aAttqAit) + (TBaqBit). (33) 

In general, the held correlation function is 




(34) 


Within the continuum models the temporal evolution of the aqi{t) is determined by dihusion 
equations. In a macroscopic description of a non-dilute solution, a dihusion current of either 
species is driven by concentration gradients of both species, so that 


JA(r, t) = -DaVca(j:, t) - DABVcs(r, t) (35) 

and 

JB(r, t) = -DBAVcA(r, t) - DbVcb{y, f), (36) 

Jj being the current of species i. In the above equations, the Di are single-species dihusion 
coefficients, while the Dij are cross-dihusion coefficients. 

Applying the continuity equation, for concentration huctuations one has 

- A^qAif) — ^ ABdqBit) (37) 
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and 


- ^ BAdqAit) — ^BdqBif)- (38) 

Here Fj = Diq^. These simultaneous equations have as solutions 

S-4(*) = [(r^- r“)exp(-r+«) + (r+ - rA)exp(-r“()] 

+ [exp(-rV) - exp(-r-«)] (39) 

and 

“«»(*) = p°T^p_ [(rB - T")exp(-r+() + (r+ - r b) exp(-r"()] 

+ [exp(-rV) - exp(-r-()] . (40) 

The predicted spectrum is 

= I+exp{—T~^t) + J_ exp(—r“t). (41) 


The spectrum of a two-macrocomponent mixture is thus a sum of two exponentials, corre¬ 
sponding to two independent relaxational modes of the system. The decay constants of the 
two modes are 


r* = ^(r4 + rB)± 


Fx-Fb 


+ ^AB^BA 


1/2 


(42) 


The in general depend on the diffusion coefficients of both solute species. 

The relaxation modes do not correspond one-to-one to solution components. Each mode 
intensity depends on the scattering powers and interactions of both species. Dehning a = 
(agyi(O)^), /9 = (aqB(O)^), 7 = ( 05 ^( 0 ) 055 ( 0 )), and A = a\a + a%l3 + 2 (Tp 1 (T 57 , the intensities 


are 

1 + = - r-)(a^a + aA(TB 7 ) + (r+ - Ta){(tI /3 + 0^057) 

+rAB(o-^7 + ctactbI^) + Tba( o-|7 -h aACBO:)] (43) 


and 

I- = [(r+ - rA)(a^a + aAfTBl) + (Ta - r-)(a| + aAa57) 

-Fab (0^7 + (^a<JbI3) - F5a(o-|7 + o-ao-bck)] . (44) 


The intensity of each mode depends on the diffusion coefficients and scattering powers of 
both species. If only one species scatters light signihcantly but both species are non-dilute. 
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the spectrum is a double exponential. If both solutes are dilute, 7 and the vanish. In 
this limit eqs. reduce correctly to the low-concentration normalized form 

a\a + a^p 

in which each exponential corresponds to the diffusion of a particular chemical species. 

D. Mutual, Self, and Probe Diffusion Coefficients 

The above results actually describe the mutual, self, and probe diffusion coefficients. To 
measure one examines a system containing a single macro component A, in which case 
eb, /S, 7 , and the F^ are zero, so F+ = F^, J+ = 1, and /_ = 0. The field correlation 
function reduces to a single exponential exp(—F^t). The continuum theory thus predicts 
that a mutual diffusion experiment on a binary system measures F^. The exponential decay 
constant Fyi/g^ may be identified with Dp or with of the microscopic theories. 

In a probe diffusion experiment, one of the components A is dilute, while the other 
component B scatters no light. Identifying A as the probe, the model requires as = 0, 
7 = 0, and F^^ = 0. The final two equalities arise because A is dilute, so almost all B 
particles are distant from any A particle and do not have their motions or positions perturbed 
by the presence of A. With these values for model parameters, the continuum treatment 
predicts F+ = F^, F“ = F^, /+ = 1, and I~ = 0, so 

t) = exp(-FAr) (46) 

Under probe conditions, according to the continuum model QELSS obtains Da of the probe. 
Da now being the diffusion coefficient of the dilute probe particles through a pseudobinary 
solution. No matter the concentration of the matrix, neither F b nor F^^ enters the spectrum. 
In favorable cases, the requirement ” ...scatters no light...” can be relaxed, as discussed below. 

In a probe experiment, fluctuations in the concentration of a concentrated matrix B may 
create currents of species A. Why don’t these currents contribute to the spectrum? Fluc¬ 
tuations aqB{0) do act on A particles, contributing a non-zero term to However, in 

the tracer limit, the currents that B induces in A are uncorrelated with the initial concen¬ 
tration fluctuations in A, so they are equally likely to enhance or diminish the fluctuations 
o-qA{t) of species A. The concentration fluctuations driven by the B-A cross-coupling are 
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not correlated with ag^(O) and do not within the continuum picture affect the spectrum 
{/(0)/(r)). 

In addition to the mutual, self, and probe diffusion coefficients, the classical literature 
includes references to a “tracer” diffusion coefficient. A tracer diffusion experiment uses 
an unlabeled species A and its labelled twin A*. The label allows one to identify the A* 
molecules, but A and A* are elsewise identical. The requirement that the two species be 
elsewise identical is more demanding than it sounds, especially for large molecules. For 
long-chain hydrocarbon polymers, even perdeuteration may lead to issues. 

In a classical tracer experiment, one creates a non-equilibrium system containing macro¬ 
scopic, countervailing gradients in the concentrations of A and A*, the gradients being so 
arranged that the total concentration ca + ca* is everywhere the same. This arrangement of 
concentrations arises in fluorescence photobleaching recovery, in which A and A* correspond 
to the bleached and unbleached molecules under consideration, with the total concentration 
Ca + Ca* being the same as the pre-bleaching concentration of unbleached and not-yet- 
bleached molecules. The flux of A* down its concentration gradient is measured, the tracer 
diffusion coefficient being 

3 a* =-DTVcA*{r,t) (47) 

As seen in the next section, the use of countervailing concentration gradients cancels any 
gradients in the non-ideal parts of the chemical potentials of A* and A. A discussion of 
the physical interpretation of the tracer diffusion coefficient resumes after a treatment of 
reference frames and irreversible thermodynamics. 


E. Reference Frames, Irreversible Thermodynamics 


This subsection considers reference frames and their implications, following closely the 
discussion of Kirkwood et al.[ll|. We discuss practical diffusion coefficients Dij and also fun¬ 
damental diffusion coefficients fltj. We also consider implications of the Onsager reciprocal 
relations, as they arise in irreversible thermodynamics. 

The practical diffusion coefficients relate the diffusion currents to the concentration gra¬ 
dients via 

g 

3i = '^DijVcj{r,t). (48) 

i=i 
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The fundamental diffusion coefficients relate the diffusion currents to the chemical po¬ 
tential gradients via 

<? 

Ji = ^ VLijViijir, t). (49) 

i=o 


The Gibbs-Duhem equation SlJ constrains the chemical potentials pj by 

djjij 


E 

i=o 


= “■ 


(50) 


Here sums pass over thermodynamic components {0, 0 denoting the solvent. 


Gibbs 


3l| emphasizes that choosing a particular component as the solvent is an arbitrary 


act. Any one component may be identified as the solvent; the other components are then 
solutes. Thermodynamic chemical components need not correspond to chemical compounds 
in a simple way, a result particularly useful in treating ionic solutions. In the following 
discussion it is convenient to use mass units, so that c is in gram-cm“^, Vi is the partial 
volume per gram of solute z, etc. 

Reference frames appear implicitly in the continuum analysis. Eq. 08] relates diffusive 
currents Jj(r,t) to local concentration gradients Vcj(r,f). Elementary physical considera¬ 
tions show that velocities (and therefore currents) only have meaning when the local zero 
of velocity—the reference frame—is specified. As a notational matter, reference frames are 
indicated by an exterior subscript R, so (J*)/? denotes the current of i as observed in the R 
reference frame. To determine a current in a frame A from the current in a frame R, one 
has 

(Ji)yl = (Jj)^ — (51) 

\ab being the velocity of frame A as measured in frame B. 

At least four frames are useful in a continuum description of diffusion. In the mass-fixed 
frame M, the center of mass of the system does not move, so 

<? 

5^(Ji)M = 0. (52) 


i=0 


Onsager[32j proposed that the fundamental diffusion coefficients are subject to symmetry 
constraints—the “Onsager reciprocal relations”—applicable in the mass-fixed frame: 




( 53 ) 
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In the solvent-fixed frame 0, the solvent is stationary: 


(Jo)o — 0. 


( 54 ) 


The solvent current in question is the bulk flow of solvent. Interparticle hydrodynamic 
interactions are hidden in the diffusion tensors. Since the fundamental hydrodynamic in¬ 
teraction tensors b and T as discussed in Section IV are derived by requiring that solvent 
flow vanishes as r ^ cxo, it is often assumed that hydrodynamic calculations are made in 
the solvent-fixed frame. Microscopic calculations of later sections refer to results in the 
experimentally-accessible volume-fixed frame. 

The relative velocity of the mass- and solvent-fixed frames is obtained by writing eq. [51] 
for component 0, and applying eq. [5H showing 


Vqm — 



Co 


(55) 


One can define a new set of fundamental diffusion coefficients (Op)o, which satisfy the 
Onsager reciprocal relations in the solvent-fixed frame if the satisfy these relations. 

From eq. [55] 

(Ji)o = (Ji)M- ^(Jo)m- (56) 

Co 

Using eq. [H] to express the (Ji)o in terms of the fundamental diffusion coefficients in 
the mass-filled frame, and applying the Gibbs-Duhem equation to eliminate reference to the 
chemical potential of the solvent, one may write 

<? 

(Ji)o = y~~^(Op)oV/ij, (57) 

i=i 


in which the (Gp)o, defined by 


(Ojj)o i^ij) M M 

Co 


^(0,o)m + ^(Ooo)m, 
Co eg 


(58) 


manifestly have the same symmetry as the 

Actual experimental data refers not to the solvent-fixed frame but to the cell-fixed frame 
c. In many systems, volume of mixing effects are small. In this case, the cell-fixed frame 
and the volume-fixed frame V, defined by 

Q 

Y,9,)vVi = 0, (59) 

j=0 


20 




are the same. The partial volume of component j is Vj. If volume-of-mixing effects are 
not small, in a classical gradient diffusion experiment the interdiffusion of the two diffusing 
components is accompanied by bulk flow and a change in the volume of the solution. 

To unite experimental data with theoretical results and with ratioanles dependent upon 
the Onsager reciprocal relations, the volume- and solvent-hxed frames must be linked. The 
relative velocity of the solvent- and volume-hxed frames is obtained from eq. [5T] by multi¬ 
plying by applying eqs. |5l] and [59l giving 


Vvo = 




i V, 


( 60 ) 


Z=1 


from which follows the relationship between the practical diffusion coefficients in the two 
frames. Using eq. 08] to replace the Jj with the Dij 




or 


k=l 

<? 


{Dij)o — {Dij)v - ^—'y^Vk{Dkj)o, i,j G (l,g) 

""0^0 fc=l 


To obtain the flj,- from the Dij, note that eq. 


can be written 

dfik. 


(J,)o= 

^ OCj 


k,l=l 


(61) 


(62) 


(63) 


Eq. 0H] then implies that the practical and fundamental diffusion coefficients obey 


(A/)o = X](^^ 


ik)0' 


k=l 


d^k 
dci ' 


(64) 


Matrix inversion techniques applied to the can be used to calculate the in terms 

of the (Zljj)o, permitting experimental tests of the Onsager reciprocal relations 11|. 

As we will see below, the reference frame correction goes away if one uses the correct 
form for the interparticle hydrodynamic interaction tensors, but the reference frame term 
connecting the solvent- and volume-hxed frames reappears in a new guise with very nearly 
the same net effect. 


F. Generalized Stokes-Einstein Equation 

We now turn to the so-called Generalized Stokes-Einstein equation, which connects the 
mutual diffusion coefficient to the self diffusion coefficient at the level of approximation of the 
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continuum models. We begin by showing that reference frame treatments serve to connect 
Dm and Dg to the flij and to the concentration dependences of the chemical potentials. 

Specializing eq. EUto a system with q = I shows {Dii)v in a binary solvent:solute system 
is the usual mutual diffusion coefficient: 


iDu)v — (-Dii)o(l — 0), 


( 66 ) 


where (j) = cihi is the volume fraction of the macrocomponent in solution. The 1 — 0 factor 
is the well-known reference frame correction to Dm- 

The self-diffusion coefficient may be measured by labeling a few solute molecules, so 
that we have a labelled species 1 and an unlabeled species 2, establishing in the system 
countervailing gradients Vci(r, f) and Vc 2 (r,f) so arranged that the total concentration 
Ci(r, f) -|-C 2 (r, t) is everywhere the same, and measuring the flux of the dilute labelled species 
1 . 

The continuum model was applied to this problem in refs. Q and y, whose treatment is 
now followed. Under these conditions, Dg may be dehned 


^ = (D,)v^V^ci(r,t). 

In terms of eq. HE] (written for two components), eq. 1661 becomes 

dci 


dt 


= ((Dn)v-(A2)y)V^ci(r,t). 


( 66 ) 


(67) 


The {Dij)v can be expressed in terms of fundamental diffusion coefficients in the solvent- 
hxed frame by applying eqs. ED and EH showing 


{Dij)v = ^ 


As seen in refs. 


y and y. 


1=1 




k=l 


djj-i 

dcj' 


( 68 ) 


the chemical potential derivatives can be written 

dni ksT 
= Oij -h I 


dc 


(69) 


'-'j '-'j 

to hrst order in Cj, 6ij being a Kronecker delta and I being an interaction integral. Similarly, 
the Qij may formally be written 




Ci 


and 


CiCj 


^ij — r 

It 


fi 


(70) 

(71) 
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fi and fij being formal dissipative factors. The concentration dependences of /i/ and flij 
ensnre that the Onsager reciprocal relations are satished, and that in a 1-component system 

'dfii 


iDii)v — (f^ii)o 


dci 


( 1 - 0 ) 


or eqnivalently that 


c — 

{D„)v ^ ^(1 - 4>)- 


/i 


On the other hand, from the above 


Dt= (O 


‘iijo 




keT 

Cl 


ksT knT 




Cl 


C2 


(O 


12j0 


and 


lim Dt = 

ci^O 


kaT 


/i 


(72) 


(73) 


(74) 


(75) 


If /i is determined by the total concentration ct = ci -1- C 2 , and if (Zlii)y and Dt are 
measnred at the same total concentration ct, eqnations [75] and [73] lead to the generalized 
Stokes-Einstein eqnation 


{Dm)v — DtCi 


dill (1-d) _ „ ( sn\ 


^BrhTT (!-'«■ 


(76) 


dc\ ksT \dc J prp 

Here H is the osmotic pressnre at constant temperatnre and total pressnre. The general¬ 
ized Stokes-Einstein eqnation reqnires that the mntnal and self diffnsion coefficients share a 
single friction factor. That is, if as explained below one writes 

i_ikBT{l - 0 ) 


D^ = [g^^\kM~ 


fu 


and 


= 


ksT 

fs ’ 


(77) 


(78) 


then according to the continnnm model the drag coefficients for Dm and Dg are eqnal, i.e., 

/m = fs- 

The continnnm model does not agree with the microscopic models given below, becanse 
the microscopic models predict that fg and /m have nneqnal concentration dependences. In 
terms of resnlts in later sections: If one expands fi = /o(l-l-aj0), the microscopic prediction 
is cis ^ aM for hard sphere snspensions. If the hard spheres gain an electric charge, adding a 
Debye potential to their interactions, microscopic models predict that the difference between 
/m and fg increases. Charging the diffnsing spheres rednces 


33( 1 Icisl towards zero, becanse 
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the spheres stay further apart, weakening the hydrodynamic forces that retard self-diffusion. 
Charging the spheres makes \aM\ larger, because the integral in eq. ll64l over the Oseen tensor 
increases. 


III. CORRELATION FUNCTION DESCRIPTIONS 


Correlation functions provide the fundamental description within statistical mechanics for 
time-dependent processes and transport coefficients. This section treats correlation function 
descriptions of diffusion, including correlation functions that describe QuasiElastic Light 


Scattering Spectroscopy (QELSS), Fluorescence Correlation Spectra (ECS) 


Q, £, 




Raster Image Correlation Spectroscopy (RICS)[3a, 1^, and Pulsed-Gradient Spin-Echo Nu¬ 
clear Magnetic Resonance (PGSE NMR38|]) The analysis here follows closely our prior 


papers, including refs. 


38 
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The reader will notice that this approach takes us from the observed relaxation functions 
through to the moments of two displacement distribution functions and then grinds to a halt. 
Calculations based on microscopic models that are powerful enough to give quantitative time 
and concentration dependences appear in a later section. 


A. Quasielastic Light Scattering Spectroscopy 


The starting point is the field correlation function, eq. [5l which gives the spectrum in 
the form an average over the positions of a particle i at time t and a particle j at time 
t + T. The average is in principle calculated as an average over a probability distribution 
function P(rj(t),rj(f -|- r)), which gives the probability of finding particles i and j (which 
may be the same particle, i = j being allowed) at the indicated locations at the times t and 
t + T, respectively. However, it is more effective to divide P(rj(f), rj(t -|- r)) into its self and 
distinct parts, and then to analyse these two parts separately. N. B.: We have switched 
back to using the to represent the locations of individual scatterers, not to label locations 
in the scattering volume. 

We begin with eq. 0 namely 


N 


g^^\q,T) = ( ^aiajexp(iq- [ri{t + t) - u(f)]) j . 


(79) 
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The terms in the sum partition into self and distinct parts, and respec¬ 

tively, with 

9W(,,t)=9<-I(,,()+J<“''(9,«). (80) 

Here 

■ [ri{t + t) -Ti 

for the self part, and 


(i)l) ( 81 ) 





exp(fq 


9'“'(9.t) 


N 


(Tj(Tjexp(fq ■ [rj(f -1- r) 



(82) 


for the distinct part. In the distinct part, the summation indices i and j are each taken 
separately from 1 to N, but the N terms in which i and j happen to be the same are omitted. 
One usefully introduces new variables: 


Xiir) = Xi(t + t) - Xi(t) 


(83) 


represents how far particle i moves parallel to the scattering vector q during time r, while 

Rij{t) = Xi{t) - Xj{t) (84) 

is the component parallel to the scattering vector q of the distance between particles i and 
j at the initial time t. 

We now consider how to evaluate or expand the two parts of the field correlation function. 
We apply the substitution 


q ■ (ri(f + t) - rjir)) ^ Vi(r) + Rijit). 


(85) 


The self part of the held correlation function may be written in terms of the new variables 


as 


N 




( 86 ) 


. 2=1 


The particles are all the same, so the sum on i can be replaced with a count N of the number 
of identical terms in the sum; the label i is now irrelevant. A Taylor series expansion followed 
by interchange of the sum and the ensemble average gives 


g 


{Is) 


(g,r) = 


i n=0 


{nX{r)y 

n\ 


(87) 
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To advance farther, one notes that the average is over the displacement distribution function 
P{X, t), this being the function that gives the probability that a particle will move through 
X during time r. Up to constants, the averages 

= j dX{iqX{TyrP{X,T), (88) 

are the moments of P{X, t). The odd moments vanish by symmetry, namely the likelihoods 
of observing displacements +Xi and —Xi are equal. P{X, r) is influenced by whatever else 
is in the system, e.g., non-scattering components, and represents an ensemble average over 
positions and momenta of those components. 

Further rearrangements lead to 

r) = Af exp + i«'‘((V(r)*) - 3(A:(r)2)2) - . (89) 

The relaxation of r) is thus determined by the even moments of P{X, r). 

A similar approach may be used to evaluate The hrst step is to introduce 

a hue-grained displacement distribution function P(X, r, {r^}). This function gives the 
probability that a particle 1 will have a displacement X during a time interval r, given that 
the coordinates {r^} of the other particles in the system at the initial time t are specihed. 
The list of M other particles includes the N — 1 other scattering particles in the system 
and the M — iV -|- 1 other non-scattering particles in the system; M = N — 1 and iV = 1 
are allowed. In studies of simple mutual diffusion, M = N] all particles are scattering 
particles. In probe diffusion experiments, N is made sufficiently small relative to the size 
of the container that interactions between the scattering particles are quite small. One can 
also study mutual diffusion of non-dilute scattering particles in the presence of a background 
non-scattering matrix macromolecule. A few experimental studies of this circumstance have 
been made; see Philliesj^, Section 11.6, for a review of these. 

The displacement distribution function P(X, r) is related to V by 

P{X,t) = j dv^ViX,T,{v^})expi-/3{WM - A)). (90) 

Here the integral is over the positions at time t of all scattering and non-scattering particles 
other than the particle of interest, while Wm is the total potential energy including the 
particle of interest and the other M particles, and A is the nonideal part of the Helmholtz 
free energy. 
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So long as the particles are all identical, t) is a snm of A^(iV —1) identical terms. It 

is convenient to make a Taylor series expansion on Xi (r) while leaving R 12 in an exponential, 




(!?Xi(t))' 


\n=0 


m 


exp(zqRi2) } ■ 


(91) 


On replacing the formal average (• • •) with the fine-grained distribntion fnnction. 




(g,r)=iV(iV-l) / 


(zqXi(T)y 


n=0 


n\ 


X j exp{-iqRi2)V{Xi,T,{r^})exp{-l3{WM - A)). (92) 

The second line of eq. |92] implicitly defines a rednced distribntion fnnction 


P 2 iXi, T, R 12 ) = J dridrupdr^dr^ ... dr^ViXi, r, {r }) exp{-(]{WM - ^)), (93) 

Here dri 2 p represents the integral over the two components of R 12 that are perpendicnlar 
to the scattering vector. We actnally need the spatial Fonrier transform 


P 2 {Xi, T,q) = j dRu exp{-zqRi 2 )P 2 {Xi, r, R 12 ), (94) 

in which Ru is the component of R 12 that lies along the scattering vector. 

The distinct part of the field correlation fnnction is finally reduced to 

= j dX,f^AXXlp(x„T,q). (95) 

Comparison of eqs. 1881 and 1951 reveals an important result. g^^^\q,T) and r) are 

obtained as averages of the same variable iqXyr) over different correlation functions, namely 
F(X, r) and P(Xi, r, g), so their information contents are not the same. We will henceforth 
drop subscripts on X when they are not significant. 

Why did we start with P{X,t) when we calculated g^^^^qp), but with F(Xi,r, {r^}) 
when we calculated g^^'^^q, r)? First, P{X, r) is obtained as an average over V{Xi, r, {r^}), 
but in that average all of the {R^} are treated the same way. When we calculate g^^'^\q, t), 
the object being averaged is Xexp{iqRi 2 )- A single coordinate F 12 is treated differently 
from all the rest, so the average must be done explicitly. Also, the count of terms in the 
intermediate step replacing Rij with F 12 was facilitated by invoking V. 
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B. Alternative Methods for Measuring Single Particle Diffusion Coefficients 


This section treats some other methods for measuring single particle diffusion coeffi¬ 
cients, namely Pulsed-Gradient Spin-Echo NMR (PGSE NMR), Fluorescence Gorrelation 
Spectroscopy (FGS), and Raster Image Gorrelation Spectroscopy (RIGS). The first of these 
measures the same as does QELSS when QELSS is used in probe diffusion mode, 

but on different time and distance scales. The other two techniques also lead back to averages 
and moments for P(X, r), but by slightly different paths. 

The physical quantity measured in PGSE NMR is quite different from the quantity mea¬ 
sured in FGS, but the mathematical form is the same, namely the PGSE NMR relaxation 
spectrum is given by 

M(2T) = M(0) exp(—iq ■ (r(t + t) — r{t))). (96) 


Here 2T is the time interval required for the formation of the spin echo, namely twice the 
time interval T between the initial 7 r /2 magnetizing RE pulse and the tt pulse that reverses 


the chirality of the spin magnetization. 


The wave vector has a new meaning, q = 7(5g, 


where 7 is the gyromagnetic ratio of the spin being observed, 6 is the duration of a gradient 
pulse, and g is the field gradient of the superposed gradient pulse. Finally, r is the time 
interval between the two gradient pulses; it is the time over which diffusion is observed. The 
notations in use the QELSS and NMR sub disciplines are not entirely compatible. For ease 
of reading we have forced eq. [96] into QELSS notation. The meaning of q in eq. [96] is very 
different from the meaning of q in the previous section, but the mathematical structure for 
the value of the PGSE NMR relaxation function is precisely the same as the mathematical 
structure was for in light scattering. The time dependence of M(2r) from the 

correlation function approach therefore looks precisely the same the expression for the time 
dependence of g^^^\q,t) using the same correlation function approach. 

Fluorescence Gorrelation Spectroscopy and Raster Image Gorrelation Spectroscopy are 
very similar in their physical description, but differ at one key point. In each method, 
one observes the motions of fluorophores, which may be free-floating molecules or may be 
fluorescent groups physically or chemically bound to the diffusing molecules of interest. 
During the experiment, a small volume of solution is illuminated with a laser beam. In 
modern instruments, the diameter w of the illuminating laser beam may be as small as a 
few hundred nm across. The fluorescent intensity is then measured at a series of times. 



and the intensity-intensity time correlation fnnction of the emitted light is measured. PCS 
and RIGS differ in that in PCS one repeatedly observes fluorescence emitted by a single 
volume of solution, while in RIGS one repeatedly observes fluorescence emitted from a series 
of neighboring volumes of solution. 

PGS and RIGS sound much like QELSS, except the re-emission process is fluorescence 
rather than quasielastic scattering. Because all phase information is lost during the flu¬ 
orescent re-emission, in PGS and RIGS one captures information on the number of fluo- 
rophores in the illuminated volume, but loses most information on the relative positions of 
the fluorophores within the volume. Qualitatively, in PGS one waits for moments when the 
fluorescent intensity is particularly bright or dim, and then waits to see how long is typically 
required for the fluorescent intensity to return to its average level. Qualitatively, in RIGS, 
one waits for moments when the fluorescent intensity from the hrst volume is particularly 
bright or dim, and then asks how rapidly the brightness or dimness spreads to neighboring 
solution volumes. 

The PGS and RIGS spectra G(r) are given by the same expression 

G{t) = f dr f dr'I{r)I{r')P{r' — r,T). (97) 

Jv Jv 

In this equation, X(r) and X(r') refer to the position-dependent intensities of the illuminating 
laser at the times t and t-|-r, respectively. Por the model calculation here, the intensities are 
approximated as being Gaussian cylinders. In PGS, the two cylinders have the same center. 
In RIGS, the two cylinders are displaced by a distance a that is for mathematical convenience 
taken to lie along the x axis. P(r' — r, r) is the displacement distribution function, expressed 
in terms of the initial and hnal positions r and r' of a diffusing fluorophore. 

Equation [HT] is a convolution integral, which can be evaluated via Pourier transform 
techniques. Por PGS, the spatial Pourier transform of the intensity 

X(r) = Joexp(—(98) 

is 

X(g) = Joexp(-gV/4) (99) 

where r is the distance from the center of the cylinder and q is the spatial Pourier transform 
vector. The spatial Pourier transform for each dimension of the displacement distribution 
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function is 


= J dx exp{iqx)P{x,T). (100) 

After a Taylor series expansion of exp{iqX), enforcing the requirement that that P{x,t) is 
symmetric in x, and extracting a leading Gaussian term exp(—g^(X^)/2) from the expansion, 
F{q, t) is found to be 

F(g, r) = exp(-g 2 ir 2 (r)/ 2 )(l + q*K^{T)/2A + q^K^{T)n2Q + (0)(g«)). (101) 

The hrst few K 2 n are 

^2(r)) = (102) 


K,{T)) = {{x{T)f)-mx{T)f)\ (103) 

and 


K,{t)) = -{{x{T)f) + 15((a;(r))2)((a;(r))^) - (104) 

In the above three equations, x(r) is the (time-dependent) displacement of the particle along 
the X axis. The brackets (...) denote an average over P{x,t). The K 2 n can all be written 
entirely in terms of the even spatial moments of P{x, r). 

Substituting into the spatial Fourier transform form of eq. [971 

G{t)=N I dq{I{q)yF{q,T), (105) 


the FCS spectrum may be written 


G{t) = 


27,11 


w'^ + K2 {t) 


1 + 


KJt) 


+ 


K.(t) 


3{w'^ + K2 {t))‘^ + K2{t))^ 

The natural variables for writing this form are the Kn{t)/w‘^^, for which 

1 K4^{t)/w^ Kq{t)/w^ 


+ ... 


(106) 


G(r) = 


27,11 


w‘^ 1 + K2 {t)/w‘^ 


1 + 


+ 


+ 


(107) 


3{1 + K2{t) / lh{l + K2 {t) / 

Many published discussions of FCS invoke the Gaussian diffusion approximation. For 
Gaussian diffusion, K 2 is the mean-square displacement, while iF 4 , Ka, and higher are all 
zero. In this case, the spectrum reduces to the form of Magde, et al.[34| 

1 


G(r) = 


2vrJ2 


1 + K2{t)/w‘^ 


( 108 ) 
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The calculation of the RIGS spectrum is extremely similar to the calculation of the FCS 
spectrum, except that in RIGS one cross-correlates the intensities observed at two times in 
two different volumes of solution. 

The RIGS spectrum may be written 

G{t) = [ dr f (ir'X(r)X(r')P(Ar, r). (109) 

Jv Jv 

Goordinates are chosen so that the origins of r and r' are at the centers of the two illuminated 
regions, in which case the particle displacement is 

Ar = r' -|- ai — r. (HO) 

with i being the unit vector along the x axis. 

For a cylindrical beam whose center is parallel to the z-axis, the intensity can be written 
in terms of spatial Fourier transforms as 

X(r) = (27r)"^y' dq^dqyIoexp{-{{q^Y + {qyY)w'^/A) exp{i{xq^ + yqy)) (111) 

and 

X(r') = (27r)“^ J dq'JqyIoexp{-{{q'J^ + {qyf)w^/A) exp{i{x'q'^ + y'q'y)). (112) 

The two components of the spatial Fourier transform vectors are (g^,, qy) and (g(,, respec¬ 
tively. 

Diffusive motions in the x and y directions are independent, so the RIGS spectrum may 
be written 

G(r) = G.(r)G,(r) (113) 

The component Gy{T) is the same as the ^/-component G{t) for FGS, namely 

/ OO 

dq,{I(q,)fF(q„T) (114) 

-OO 

while the x-component is 

/ OO 

dg(X(g))^X(g,r) cos(ga). (115) 

-OO 

The integral for the x-component gives 

GJr) = N (exp (^- (2 (A',(r) + 
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X 1 + 


K4{t) - 6a^ {K2{t) + w"^) + 3 {K2{t) + 


24 (ir2(r)+w2)^ 

Keir) - 15a^ (i^^2(r) + w^) + A5a^ {K2 {t) + w^f - 15 {K2{t) + 

720 {K2iT)+w^f 

This form is not quite the same as that seen in Digman, et al. 36 


( 116 ) 


37| because: Digman, et 
ah, only allowed for Gaussian Diffusion; the displacement a between the two laser beams 
is here taken to be a continuous variable; and the time displacement t is here used as an 
independent variable relative to a. The observed spectrum will be a product of eq. 11161 and 
the FCS relaxation forms for the y and z directions. 


C. Partial Solutions for g^^\q,T) and g^^^\q,T) 


This section treats partial solutions in which g^^\q,T) is calculated up to some point in 
terms of molecular parameters. We start with the simplest approaches. 

The hrst case refers to systems in which the N scattering particles are highly dilute. 
There may also be M — iV + 1 matrix particles, which may not be dilute. If the scattering 
particles are adequately dilute or elsewise noninteracting, then vanishes. The 

disappearance is seen in eq. HH in which the factor exp(*gi?i 2 ) refers to a distinct pair of 
scattering particles. If particles 1 and 2 essentially never interact with each other, then 
all phases of the complex exponential are for all practical purposes equally likely, in which 
case the complex exponential averaged over all pairs of scattering particles averages to zero, 
taking the integral with it. In this case, then, g^^\q, t) is very nearly determined by g^^^^q, t), 
eq. |89l 

At the bottom end for simplicity, for dilute probes particles in a simple, low-viscosity 
liquid, P{X,t) is a Gaussian in X, in which case in eq. IHUthe terms in g®, and higher 
all vanish. The held correlation function becomes 


9'“>(9.r) = Afexp(-,2((A'(T))2)/2). 


In such a system, Doob’s Theorems 


42[ | guarantee that 


((A'(r))^) = 2Dt 


(117) 


(118) 
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so that the spectrum decays as a single exponential in time. Further details are found in 
Berne and Pecora Q, Chapter 5. 

A word of caution: Berne and Pecora were writing for an audience with a solid knowledge 
of the physics and theoretical chemistry involved. Their Chapter 5 gives a treatment that 
only applies to dilute particles in simple solvents. They reasonably expected that their 
audience would recognize this. Berne and Pecora also treat results relevant to complex 
fluids in parts of their Chapters 10-12. 

A claim that P{X, r) is a Gaussian is often said to arise from the Central Limit Theorem, 
which gives the consequences if r(f) is composed of a large number of small, independent 
steps. Instead, r{t) is all or part of a single step. Under the conditions in which the Central 
Limit Theorem is valid, so that P{X,t) is a Gaussian in X, eq. 11181 is equally sure to be 
valid. Therefore, whenever the Gentral Limit Theorem is applicable, g^^^\q,T) is sure to be 
a single pure exponential in r as well as in 

For the last four decades, there has been interest in studying optical probe diffusion 
in complex fluids, fluids in which there are relaxations on the time and distance scales 
accessible to experimental study with QELSS. However, if the fluid has relaxations on the 
time scale being studied, or even longer time scales, then r(f) is emphatically not composed 
of large numbers of independent steps. The Gentral Limit Theorem is totally irrelevant 
to measurements being made on time scales on which the fluid has relaxations. In this 
interesting case, P{X, t) may well not be a Gaussian, so (in the absence of an independent 
direct measurement showing that P{X,t) is a Gaussian, for reasons other than being a 
Gentral Limit Theorem Gaussian) eq. 11171 can not be invoked. 

Indeed, there is extensive experimental evidence showing that P{X, t) is not Gaussian in 


complex and glassy fluids on time scales on which there are ongoing relaxations. 


m, 


43N45|: 


First, particle tracking measurements can in some systems determine P(X,t) directly, at 
least on longer time scales. Fifteen years ago, Apgar, et al. [4^ and Tseng, et al.[47| measured 
P{X,t) directly for probes in polymer solutions, unambiguously Ending that P{X,t) had 
a non-Gaussian form. Wang, et ah 48|, l49|, and Guan, et ah ^ measured more precisely 
the nature of the non-Gaussian deviations in P{X,t). For small X{t), P{X,t) is close to 
a Gaussian, but at larger | X(t) |, P{X,t) only decreases exponentially with increasing 
\ X{t) I. Ghaudhuri, et ah 29j report similar Tat tails’ for glassy systems. While some of the 
experimental systems were quite complex, Guan, et al.’s system^ was physically simple. 
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namely it was a suspension of smaller colloidal hard spheres diffusing through a nondilute 
suspension of larger hard spheres. 

Second, if P{X,t) is a Central-Limit-Theorem Gaussian, log( 5 '*^^^^(g, t)) must be linear 
in and in t. For dilute probe particles in complex fluids such as polymer solutions, 
neither of these assumptions is true in general. For example, for polystyrene latex probes in 
hvdroxvpropvlcellulose:water 5]]- 53|. log(gd*)(g,t)) is a sum of two or three stretched and 


unstretched exponentials. Furthermore, for some modes but not others, the mean relaxation 
rate is not linear in 

Third, supporting evidence that P{X, t) is not a Central-Limit-Theorem Gaussian is 
provided by FCS measurements. On one hand, transverse diffusion of macromolecules in 


cell membranes was studied by Wawrezinieck, et al.[^ They found in their systems that 
the relaxation time td was linear in but did not trend to zero as w —)■ 0, so the diffusion 
coefficient inferred from td depended on w. Their measurements show directly that P{X,t) 
was not a Gaussian in X{t). Furthermore, as shown by Schwille, et al.j^|, in complex fluids 
the mean-square displacement inferred from eq. 11081 does not always increase linearly in t. 

As a practical experimental aside, if probe diffusion is nouGaussian, eqs. 18^ ITOTt and lllOl 
are still valid. By measuring at multiple scattering vectors g, or by measuring G{t) 

for a series of spot diameters w, K 2 {t), and higher-order terms may be accessible as 

functions of time. 


IV. GENERALIZED LANGEVIN EQUATION 

A. Introduction 

This Section considers how the concentration dependence of diffusion coefficients can be 
obtained from a Generalized Langevin equation approach. The calculations here represent 
extensions of the correlation function descriptions of the previous section. Because we insert 
mechanical models for particle motion and concrete forms for the direct and hydrodynamic 
interactions between the particles, we can make quantitative predictions for the dependences 
of Dra, Ds and Dp on solute concentration and other solute properties. 

In Section IIVBI I show how various diffusion coefficients are related to the hrst cumulant 
of a relaxation spectrum. In order to interpret the integrals in Section IIVBI a model for 


34 
















particle motion that properly represents the forces between the particles must be supplied. 
This representation is the generalized Langevin equation seen in Section IIV Cl Section IIVDI 
treats direct interactions between the diffusing particles, as treated with their interparticle 
potential energies. Section IIV El treats hydrodynamic interactions between particles, as 
represented by hydrodynamic interaction tensors. Claims that there is screening of the 
hydrodynamic interactions between macromolecules in solution are considered and refuted. 
Sections IIV FI and IIV Gl evaluate our expressions for Dm, including (Section IIV Gjl simple 
effects of direct and hydrodynamic interactions and (Section lIVHj) dynamic friction terms. 
Section liV Il considers the long-range behavior of the Oseen tensor, in particular the difference 
between inhnite and closed-system behaviors, leading to a microscopic replacement for the 
reference frame discussion. Section liV Jl examines the wavevector dependence of Dm- Section 
lIVKI evaluates the self-diffusion coefficient Dg and the probe diffusion coefficient Dp. 

B. Diffusion Coefficients from Cumulants 

This section presents the general formulations that lead to extracting the mutual, self, 
and probe diffusion coefficients from the first cumulants of and various 

numerical transforms, as obtained using QELSS, FPR, PGSE NMR, FGS, RIGS, and related 
techniques. It should be emphasized: Gumulants are obtained from the short-time limits of 
various time derivatives. However, cumulants represent weighted averages over relaxations 
that decay on all time scales, short and long. Suggestions that cumulant series cannot rep¬ 
resent multimodal relaxations are incorrect. However, in some cases alternative expansions 
may be more useful. 

Our starting point is eq. [Tl]for the mutual diffusion coefficient, which may be rewritten 
as 

Dmq^ = (119) 

r-i-o dr 

For simplicity, in this section we take at = 1. From eqs. lH [121 and [14] the mutual diffusion 
coefficient for identical particles is 

1 d / ^ \ 

= ~ q^gm{q,0) ^ \ ~ ^ 

N is again the number of scattering particles in the system. Applying the identity rj{t + T) = 
rj(t) -f- ds Vj(s), Vj(s) being the velocity of the particle at time s, followed by a Taylor 
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series expansion in q ■ Vj , one has 


Dm = 


d 


N 


ds Vj(s)]' 


0) r^o Ot 


*J=i 


n=0 


n\ 


( 121 ) 


where Yij = Vj{t) — Vi{t). Taking the derivative, the leading terms in q are 


N 


Dm. 


g^5fd)(g, 0) T-i^o 


lim ( y exp[iq ■ r^ht)] (iq ■ v, (t + 

T—^0 \ f ^ 


T] 




rt+T 


-q 


2 . 


( 122 ) 


ds Vj{s)\j(t + r) + ...^ ^ 

N is again the order of the £t. The two velocities in the integral refer to the same particle. 
The above steps nse only standard methods of freshman calcnlns: differentiation, integration, 
and expansion. No appeals to models of particle motion, time scales, or statistico-mechanical 
stationarity were invoked. 

The process of linking Dm to particle motions is readily dnplicated for Dg and for Dp, 
which are each the hrst cnmnlant of r), albeit in different systems. One may envision 

measuring Dg by tracking an isolated particle in a uniform solution. [A QELSS measurement 
in homodyne mode requires that at least two particles be present; elsewise S{q,t) is a 
constant.] For an experiment that tracks a single particle 


and 


g'’(q,T) = /^exp(!q- |rj(«) - r,)! + t)])\ 


D,r/ = - lini A In r)). 


(123) 


(124) 


The calculation of Dp is notationally more complicated, because the probe and matrix 
particles may be entirely different in their natures, but eq. 11241 applies to both calculations. 
That calculation is postponed until Section lIVKI 

C. Generalized Langevin Equation 


To make further progress, an adequate description of particle motion in solution is needed. 
Here particle motion will be characterized with a generalized Langevin equation. The orig¬ 
inal Langevin equation}^ described an isolated particle in an external potential 

rdv{t) 


M- 


dt 


= -Mt) + FB{t)-VW{r,t), 


(125) 
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where M and f are the particle mass and drag coefficient, and W(r, t) is the potential energy 
of the Brownian particle in the external held. Vhh may depend slowly on position, bnt only 
on length and time scales far longer than those over which varies. 

Farther interpretation of eq. 11251 reqnires a discnssion of the signihcant time scales in 
the problem. If one applies a force to a particle, there is some short time ~ th before the 
snrronnding solvent molecnles reach their steady-state behavior. Only for times f > r// is 
Stokes’ Law behavior or anything similar expected. A second scale tb = rn/f describes 
the time reqnired for inertial relaxation of the Brownian particle. Rice and Gray ^ show 
Tb > Th- 

Over times t ^ tb, which are the only times nsnally accessible to qnasi-elastic light 
scattering, 0. For t S> tb , the particle velocity v in the Langevin eqnation may be 

divided as 

w = wb + ^d- (126) 

The drift velocity is defined as 

VD = -r^VWiv)- (127) 

V£) is the sedimentation velocity fonnd in the Smolnchowski sedimentation-diffnsion eqna¬ 
tion. 

The Brownian velocity describes the solvent-driven motions of the particle over times 
t > Th- We are considering real diffnsing particles, not solntions to stochastic eqnations, so 

is a continnons variable with well-behaved derivatives and integrability. arises from 
stress flnctnations in the solvent and at the solvent-particle interface, as described by the 
random force FB(f). is nsnally assnmed to have a very short correlation time ~ tb'- 


(Vb(«) ■ Vb(0)) = 0, 


(128) 


There are some technical complications related to v(t) being unsteady, namely that Stokes’ 
Law must be replaced by the Boussinesq equation. For a discussion of this point, see Chow 
and Hermans M |. 

Unlike v^, may have a slow secular time dependence, but is virtually constant over 
times t ^ Tb- Random motions are generally assumed to be decoupled from secular drift 
velocities v^j, so that the Kubo relation 

^ ds vb(s)vb( 0)\ = DI, f > Tb, (129) 
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is assumed to continue to apply if the particle is in a near-constant external potential. Here 
I is the identity tensor. 

The Langevin equation was originally applied to diffusion in an external field, such as 
the helds encountered in a centrifuge or electrophoresis apparatus. In these systems, the 
Langevin approach and the Smoluchowski diffusion-sediment at ion equation are equivalent, 
as seen in standard textbooks[^. However, to apply the Langevin equation to compute the 


dynamics of mesoscopic particles in solution, one must make fundamental reinterpretations 
of several terms of the equation: 

First, the external potential W{r) must be replaced by an iV-particle interparticle po¬ 
tential W]sf{r^), where = (ri, r 2 ,..., r^v). This reinterpretation, which appears in the 


literature as far back as the Kirkwood-Risemann [58| papers on polymer dynamics, is not 
physically trivial. Unlike an external force VW (r), the intermacromolecular force ViWN{r^) 
on particle i is strongly and rapidly time-dependent. Variations in the interparticle forces 
are correlated with Brownian displacements of individual particles, because Brownian dis¬ 
placements Atb = f ds vb(s) in the positions of individual particles are large enough that 
W{y^ + Ay^) 7 ^ WN{r^). Furthermore, Brownian-motion-induced changes in IFAr(R'^) 
occur on the same time scale as Brownian motion. 

Second, the position-independent drag coefficient / must be replaced with an iV-particle 
mobility tensor /i^ = pjj(r^). / may be concentration-dependent because the particle of 
interest interacts hydrodynamically with other particles in the system, but in eq. 11251 only 
an ensemble-average (“mean-held”) value of / was used. In contrast to /, which depends 
only on macroscopic concentration variables, fiij depends on the current (and, in visco-elastic 
solvents not treated here, previous) positions of the other particles in solution. In a 
system of interacting Brownian particles. Brownian motions and drift velocities are partially 
correlated, leading to complications with eq. 11291 as discussed below. 

Having made this fundamental reinterpretation of the Langevin equation, for a system 
of N interacting Brownian particles, eq. 11261 becomes the Generalized Langevin Equation 


Vj = VBi + ^Di- 


(130) 


In this equation, each particle i has its own velocity Vj and Brownian and driven velocity 
components and This equation looks exactly like eq. 11261 but the above fundamental 
reinterpretations apply. We now consider what forces drive and \Di, and how they are 
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correlated. 


D. Interparticle Potential Energies 


V£)j is in part determined by the interparticle potential energy, and in part determined by 
the iV-particle mobility tensor discussed in the next section. The total potential energy W 
of the system is typically written as the sum of pair potentials V. In real physical systems, 
three-body potentials that cannot be written as sums of pair potentials are undoubtedly 
important. For simple hard spheres of radius a the pair potential energy is 


I/(r) = 0, |r| > 2a 
I/(r) = cx), \r\ < 2a. 


(131) 

(132) 


Formally, V (r) for a hard sphere lacks a derivative at r = a, so the force between two hard 
spheres is either zero or undehned, threatening mathematical complications. Physically, the 
force between a pair of particles is well-behaved everywhere; these complications should not 
arise in real systems. Final expressions for Dm involve integrals over radial distribution 
functions g{r), forces not appearing explicitly in the equations. g{r) for a spherical macro¬ 
molecule with no long-range interactions is not quite identical to g{r) for hard spheres, 
because a real macromolecule is slightly compressible and delivers a well-dehned force dur¬ 
ing a close encounter. The hard-sphere form for g{r) should still be a good approximant to 
a real g{r). 


For charged hard spheres (often used l60|] in theoretical models to approximate the 


behavior of small micelles) in solutions that also contain added salt, a screened Coulomb 
potential 


l/(r) = 


Q‘e 


2 p—K.{r—2a) 


(133) 


47rer(l -|- naY 

is added to eq. 11311 This equation gives the interaction between a single dielectric sphere 
and a point charge. Here Q is the sphere charge, a is the sphere radius, k is the Debye 
screening length, and e is the ratio of the dielectric constants of the solvent and the sphere 
interior. Note that the potential energy between a dielectric sphere and a point ion is not 
the same as the potential energy between two charged dielectric spheres; see Kirkwood 


and Schumaker 


6ll | and this author 62|. For objects the size of micelles induced-dipole 


potential energies may be as important as the potential energy described by eq. 11331 [62 1 
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Some treatments of micelles also incorporate an intermicellar van der Waals potential, e. g., 


refs. 
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E. Hydrodynamic Interactions; Hydrodynamic Screening 


In addition to the direct intermacromolecular interactions described by macro¬ 

molecules in real solutions encounter two sorts of hydrodynamic interactions. First, when 
a particle is driven through solution by an outside force, the particle sets up a wake in the 
liquid around it. The wake drags along nearby particles, so applying a direct force to one 
particle causes nearby particles to move. Scattering of the wake by other particles in solu¬ 
tion leads to higher-order many-particle hydrodynamic interactions. Second, because there 
are hydrodynamic interactions, the Brownian forces and velocities of nearby Brownian 
particles are strongly correlated, so that for i ^ j 


( j ds VBj(s)vBj(0)\ = Dij 7^ 0 , f > Tb- 


(134) 


Dij is a two-particle diffusion tensor. In general, Dij = Dij{r^) depends on the relative 
positions of all particles in the system. Our eq. 11221 only uses the i = j terms. 

From the Kubo relation (eq. 113411 one generally infers 

^ exp[iq ■ rjj(f)] / ds q-VBj{s)c[ -VBj{t + t) ) = 




N 


(135) 


^ exp[lq ■ ry(l)]q-■ ^ ), 

Eq. 11351 embodies three implicit assumptions. First, exp[iq ■ rij(t)] is taken to be essentially 
constant during the short interval within which J ds "VBj{s)\Bj{t + t) is nonzero. Second, 
the subtle correlations between and \Bj are assumed to decay between t and t + r, so that 
rj{t) and "VBj{t -|- r) are not correlated. Third, in evaluating this equation, r is sufficiently 
short that Djj can be evaluated using the initial particle positions. 


As shown in ref. 


63l |. while rj{t) and Vj(t) are required by statistical mechanics to lack 


equal-time correlations, one may not assume that -|- s) and -|- r), the latter 


depending on the particle positions at time t + r, are similarly uncorrelatec 
account of correlations between Brownian and later driven velocities leads 
expressions for D^- 


. Failure to take 


63l | to erroneous 
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Finally, 


N 


^ exp[iq ■ rij(t)]iq ■ + r) \ = 0, r'> tb, (136) 

WJ=i / 

because between t and t + r the Brownian velocity will have thermalized. At the time t + r, 
the Brownian velocity "Vsjit + t) will have lost all correlation with rij{t). 

The direct velocity of particle i, as modihed by hydrodynamic interactions, is 


N N 

VBi = ^ ^ p-p ■ Fj = — ^ ^ pp ■ VjhFjv(r ) 
i=i i=i 


(137) 


In eq. 11371 Fj = — VjTFAr(r^) is the direct force on particle j, while pij is a hydrodynamic 
mobility tensor the allows a force on particle j to create a motion of particle i. If r is 
small, is nearly constant over {t,t + r), so ds V£)j(s)v£)j(t + r) is linear in r. The 
drift velocity depends on the time-dependent relative positions of all particles, so ^rBi(t + r) 
depends on r. 

For hard spheres in solution, there are explicit forms for the p^j. From the work of 
Kynch64], Batchelor 65], this author 66], and Mazur and van Saarloosj^, p^ can be ex¬ 


panded as 




\ m,m^iorl-,m,neqi,l 


^iml H“ • • • 


for the self terms and 


Pij — 


J. [ Tjj -|- Timj + 




for the distinct terms. 


The leading terms of the b and T tensors arej^ 

, 15 / ay_ 

Oil = —r ~ 

4 \riij 


(138) 


(139) 


(140) 




75a^ 


1 ( 3.^2 „ 2„2 
im il ml 


{[1 - 3(ri,n • - 3(f^; • Ylif] 


i ' ^im ) }f im? 


(141) 


T-. = --^11 f-T'-l 
^ ' ij 


( 142 ) 
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TjmZ = 


15 

Q ^2 j,2 

^ ‘ im‘ml 


[I 3(r2m ■ ^m/) l^im^mh 


(143) 


where a is the sphere radius, rij is the scalar distance between particles i and j, rr denotes 
an outer (dyadic) product, and only the lowest order term (in -) of each tensor is shown. 


See Mazur and van Saarloos 


67l | for the higher-order terms. 


hij and Tij represent interactions between pairs of interacting spheres. T^- describes 
the displacement of particle i due to a force applied to particle j, while describes the 
retardation of a moving particle i due to the scattering by particle j of the wake set up by 
i. Timi and himi describe interactions between trios of interacting spheres. Tjmz describes 
the displacement of particle i by a hydrodynamic wake set up by particle I, the wake being 
scattered by an intermediate particle m before reaching i. hi^i describes the retardation of a 
moving particle i due to the scattering, hrst by m and then by of the wake set up by i. The 
effect of pair interactions on Dm was treated by Batchelor 65(] and a host of others; effects 
of three-body interactions were treated by this author 6^ and virtually simultaneously by 


Beenakker and Mazur 


68( |. An approach similar to that used to generate eqs. 1 140111 43 1 allows 


calculation of the hydrodynamic interactions between 2, 3, or more random-coil polymers, 
as well as the concentration dependences of various polymer transport coefficients. 69]. 

In terms of the above, the diffusion tensor is 


^ij J • 


(144) 


It is sometimes proposed [70| that hydrodynamic interactions in many-body systems 
should be screened, i. e., T ~ exp(—Kr). The term screening is meant to suggest 
an analogy with electrolyte solutions, in which the Coulomb potential has the form seen in 
eq. 11331 when a background electrolyte is present. Snook et ah [70t| discuss how screened 
hydrodynamic interactions would modify Dm- The basis of the assertion that there is hy¬ 
drodynamic screening is the correct observation that the hydrodynamic interaction and the 
Coulomb interaction are both 1/r interactions, the correct observation that the Coulomb 
interaction is screened in ionic solutions, and the incorrect conclusion that by analogy the 
hydrodynamic interaction must also therefore be screened. The analogy is immediately 
flawed in that the Coulomb interaction is a 1/r potential, while the Oseen interaction is 
a 1/r force, the longest-range force in nature. The falsity of the analogy is shown by the 
existence of gravity, as quantitatively described by Isaac Newton. Gravitating bodies have a 
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bare 1/r interaction, but a man y-p article gravitational system very certainly does not show 


screening of gravitational forces 


nl- 


Rigorous analyses of many-body hydrodynamics show conclusively that hyc 
teractions between Brownian particles are not screened. Beenakker and Mazur 


ro dynamic in- 


72| resummed 


the complete many-body hydrodynamic interactions tensors, showing that intervening par¬ 
ticles weaken the fundamental hydrodynamic interaction Tij but do not alter its range from 
1/r. The resummation included all ring diagrams, ring diagrams being the longest-range 
many-particle diagrams. While Beenakker and Mazur did not include all diagrams in their 
resummation, rigorous mathematical analogy with plasma theory shows that the more com¬ 
plex diagrams not included in their resummation are shorter-ranged than the ring diagrams, 
and cannot reduce the fundamental range of T below r~^ 

For Brownian particles, deGennes 


73 


of hydrodynamic screening. Reference 


is sometimes cited as demonstrating the importance 


73| in fact only cites Freed and Edwards |74l| . ra ther 


than presenting a derivation, and must be read against statements earlier in ref. 


73 


that 


translation invariance assures that there is no hydrodynamic screening at long distances in 
macromolecule solutions. Ref. 0 has since been considerably rehned by Freed and Perico 75], 
who conclude that hydrodynamic interactions are modihed by polymers that are free to move 


through solution, but tha 


Altenberger, et ah 


76 


hydrodynamic interactions are not screened. 


77| present an extended discussion of hydrodynamic screening 


m 


different systems. They show that hydrodynamic interactions in solution are not screened, 
at least on time scales sufficiently long that inertial effects may be neglected. The absence 
of screening can be understood in terms of momentum conservation. In a sand bed, in 
which particles are held hxed, T is known to be screened. Friction irreversibly transfers 
momentum out of the fluid flow into the sand matrix, the irreversible loss of momentum 
being directly responsible for hydrodynamic screening. In a solution, friction can still transfer 
momentum from the fluid into suspended macromolecules, thereby slowing the fluid and 
accelerating the macromolecules. However, in an equilibrium system, the macromolecules’ 
velocity distribution is independent of time. Any momentum transferred from the fluid 
into solute macromolecules must on the average rapidly return to the fluid, so interactions 
between the solvent and dissolved macromolecules do not transfer momentum irreversibly 
from the fluid to the macromolecules. On the time scales that are accessible with QELSS 
or related techniques, hydrodynamic screening therefore does not exist in macromolecule 
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solutions. Only if one omits the return of momentum from the particles to the fluid can one 

n 

obtain screening for hydrodynamic interactions in solution. Altenberger, et ah [77[ present a 
detailed analysis of all claims, up to their time of publication, that there is hydrodynamic 
screening, showing why each of these claims is incorrect. 


F. Application of the Model 


Eq. 11221 dehnes as the limit of a derivative as r —)■ 0. This section considers that 
derivative as applied to a system that follows the Generalized Langevin equation, which for 
a suspension of interacting macromolecules becomes 

N 

Vj = + /in ■ (-VifEjv) + ^ /ip(-VjTEAr) (145) 

It must be recognized that digital correlators only measure for r S> tb- The 

experimentally accessible short-time limit of only reflects system behavior occurring well 
after the Brownian velocity components, other than the long-time tail, have been thermal- 
ized. The limit r —)■ 0 is thus actually a limit r —)■ e for some e S> th,tb- Combining eqs. 
1122111301113411135111361 and 11371 and eliminating terms that are linear or higher in r one 
hnds as an intermediate form 


Dm 


l>t+T 


g2^d)(g,0) 


( (^exp( 2 q ■ rij(t))(-iq ■ - zq ■ voi 




+ J ds exp[iq ■ rjj(t)](q ■ [vBjis)vBjit r) -h ^Djis)vBjit + r) 
+ VBj(s)vB,j(t -h r) -h vj:)j(s)vB,j(t -h r)] ■ q) ) ). 


(146) 


The hrst and last terms of the series of terms in vbj and vb vanish, the hrst because 
between t and t + r VBj(t + r) has thermalized and the last because it is linear in r and 
vanishes as r ^ e. The fourth and hfth terms, the terms in vn,VRj are the dynamic friction 
contribution to Dm, (the name was hrst used by Chandrasekhar 78(] in connection with a 
similar problem in stellar dynamics), whose individual terms are 

q ■ AD^ ■ q = ds exp[iq ■ rj^jq • [vBj(s)vDj(t + r) -F VBj(s)vBj(t + r)] ■ q^ . (147) 
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and whose collective effect is 


N 


q ■ ADm ■ q = q ■ ■ q 


(148) 


ADm describes the correlations between the Brownian displacements of each particle 
and the snbseqnent direct forces on it. AD^ is non-vanishing becanse, on any time scale 
snfhciently long that S'(g,r) is not a constant, the particles have moved from their initial 
positions. The expression for rednces to 

1 


Dm = 




^^^exp(iq-rij(t))(-iq-Voi-q-Ai-q)^ - q ■ ■ q^ . (149) 


An integration by parts displaces the points of application of the Vj operators implicit in 
V£)i; the Vj were shown explicitly in eq. 11451 One finds 

1 


g2^(l)(g,0) 


I keT I ^ iqVi ; [/i^z exp( 2 q ■ r^)] - ^ q ■ ^^jj ■ qexp[iq ■ Vij] ) - q ■ ADm ■ q 

\ \ i,j / 


(150) 

Eq. 11491 represents most clearly that Dm arises from two sorts of motion, namely self- 
diffnsion of individnal solnte molecnles (terms in and ADm) plus driven motion dne to 
interparticle interactions (the term in v^ij). Drift velocities are typically mnch smaller 
than Brownian velocities v^j snggesting (incorrectly!) that the self-diffnsion terms shonld 
dominate driven motion in eq. 11501 However, drift velocities persist over long times, while 
Brownian velocities at a series of times add incoherently to the total particle motion. What 
is the persistence time of the drift velocity associated with a concentration ffnctnation aq{t)l 
So long as aq{t) persists, the macromolecular concentration will be nonnniform. A macro- 
mo lecnle in the concentration gradient corresponding to aq{t) will on the average experience 
a net force that persists as long as the concentration gradient does. On the time scales 
observed by QELSS, the drift contribntion to Dm can be several times the self-diffnsive 
contribntion. As an experimental demonstration, note Doherty and Benedek’s demonstra¬ 
tion that removing the salt from a sernm albnmin solntion can increase Dm as measnred by 


QELSS by severalfold over its high-salt valne 
Eq. 11501 simplifies. Ref. 

1 


63 


79|. 


demonstrates 


Dm = 




exp[iq ■ r^jiqVz : [Dji] + ^^(q ■ Dij ■ q) exp[zq ■ Vij] ) -h q ■ ADm ' q j 
i,j,l / / 


( 151 ) 
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As emphasized by Felderhof 80|, V; ■ Dji 7 ^ 0 if Dji is taken beyond Eq. 11511 is 

totally consistent with the Smoluchowski equation if one interprets D of the Smoluchowski 
equation, not as the Kubo-form coefficient of eq. 11341 but as a dressed diffusion coefficient 
incorporating the dynamic friction correction 63|] q- ADm ■ q. Hess and Klein’s 811] excellent 
review presents a different treatment of this same question. 

The evaluation of eq. 11511 naturally separated into two parts. In the next section, we 
evaluate the explicit ensemble average (■■■). The following section treats the dynamic 
friction term q ■ AD^ ■ q- 


G. Evaluation of Dm for Hard Spheres 


For a solution of hard spheres, the most extensive evaluation of eg. 11511 is that of Carter 


and Phillies 8 ^; calculations of Dm are also reported by Felderhof 


and Beenakker and 


Mazur [68(1 . Eq. 11511 can be written as a cluster expansion, i. e. as a sum of averages over 
distribution functions. The n-particle distribution function is 




-N-j-n 


d{N — n} exp[—/ 5 (hF 7 v — A)], 


( 162 ) 


where V is the system volume, A is the normalizing coefficient V~^ J d{N} exp[—(]{W^ — 
A)] = 1, and d{N — n} denotes an integral over macroparticles N — n+l,N — n + 2,.--,N. 
The being concentration-dependent, have (except in electrolyte solutions) pseudovirial 
expansions 

+ c^g^^’‘^\ri2) + ■■■ (153) 

g^^\ri 2 , ris) = ris) + cg^^’^'>{r 12 , ris) + ••• (154) 

In these expansions the are concentration-independent. 

The terms of eq. I151l that contribute to Dm through second order in c may be rearranged 
as 


Dm = 


Do 


\Nq^ + Ia + Ib Ig] 


in which 


g2^(b(g,0) 

Ia = cN j dr 5 '^^Hr)q ■ bi 2 (r) ■ q 
Ib = cN j dr 5 f*'^^(r)q ■ Ti 2 (r) ■ qexp(—iq ■ r) 


(155) 


(156) 


(157) 
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Ic = -cN j dr g^‘^\r)[exp{-iq ■ r)iqV : [bi 2 (r) + Ti 2 (r)] + iqV ; Ti 2 (r)] (158) 

Id = c^N J (ir(is s)q ■ bi 23 (r, s) ■ q (159) 

Ie = c^N j (irrfs s)q • Ti 23 (r, s) ■ qexp(—-iq ■ s) (160) 

If = —c^N j (ircis s) exp(—iq ■ s)iqV : [bi 2 (r) + Ti 2 (r)] (161) 

Ig = -c^N j drds s)[l+exp(-zq-r)+exp(-iq-s)][iqVi : bi23(r, s)+iqV3 : Ti 23 (r,s)] 

(162) 

where r = ri 2 and s = ri 3 . Replacing the with the ^("’0 of eqs. 11531 and 11541 transforms 
the li into a psendovirial expansion for Dm- 

In the limit of low g, most of these integrals were evalnated by Carter 82(] using pre¬ 
dominantly analytic means or by Beenakker and Mazur 68(] by Monte Carlo integration; 


Chere are small differences 
68l. l83| omit the V ■ /i terms, 


Felderhof 83j] obtained many of the corrections analytical 
between different calculations of Dm because some authors 
while others treat reference frame corrections in different ways. Table I, taken from Carter 
et ah 82(] gives the contributions of the Jj to the first and second order concentration cor¬ 


rections to Dm- In Table I, the fundamental concentration unit is the solute volume fraction 
(j) = 4:7ia^N/3- 

Table One 

Integrals for the concentration dependence of Dm = Do[l + Kicj) -\- /c20^]/5'^^^(g, 0). 


Equation 

Ki 

h 

Ia 

-1.734 

-0.927 

Ib 

-5.707 

13.574 

Ic 

-1.457 

1.049 

Id 

0 

1.80 

Ie 

0 

6.69 

Ip 

0 

-4.1348 

Ig 

0 

4.120 

totals 

-8.898 

22.17 


Details of the integrations appear in Carter et al. 
worthy: 


. Several technical issues are note- 
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1) In eq. 11571 the convergence of over the (a/r)^ component of Tjj is delicate. The 
Oseen tensor “(I + ff) is only an approximation to the long-range part of the exact T'^, 
inaccurate at extremely large distances because it does not enforce the physical requirement 
that the total volume flux of solvent and particles, across any surface that divides the 
container in twain, must vanish. Correction of the long-range behavior corresponds to 
the reference frame correction [h|, which physically requires ~ for a system with negligible 
volume of mixing - that the volume flux of solute into a volume of space must cancel the 
simultaneous volume flux of solvent out of the same volume.The next subsection expands 


on this point. As shown in refs. 


66 


84| and the next subsection, for a small hard solute the 


reference frame correction may be stated 


dr ([q ■ T’(r) ■ q] exp( 2 q ■ r) -|- (j)iH{q)) = 0 


(163) 


which lets one write 


q 


dr 5 f^^’°^(r)T’(r) = q^ 


dr 


(^( 2 , 0 ) ^ _ 0^^2 


(164) 


where (pi is the fraction of the system’s volume occupied by a single particle, Thiru) is the 
true long-range part of the hydrodynamic interaction tensor for a closed container, H{q) is 
the Fourier transform of the hydrodynamic shape of the solute particle, and T’(r) has been 
approximated by T(r) over the short range within which g{r) — 1 ^ 0. Eq. 11631 is basically 


physically equivalent to Oono and Baldwin’s 


method for removing the divergence from 


the perturbation series for of a solution of random-coil polymers. 

2) A term of eq. 11571 is 

I r . /AN 3 


Ib 2 = -c dr g^^’°\r)[q^ - 3(q ■ f)] exp(- 2 q ■ r) 


(165) 


where the limits on |r| come from If the exp(iq ■ r) were absent, the angular integral 

over [g^ —3(q-f)] would vanish, while the J°° r^dr (a/r)^ would diverge; the integral without 
the exp(iq ■ r) is improper and has no meaningful value. On retaining the exp(iq ■ r) term, 
doing the integrals, and taking the small-g limit at the end, one finds Ib 2 = 4:7ca^cq‘^/3. 

3) To evaluate Id-■■ da analytically. Carter et ah 82|] used a spherical harmonic expan¬ 
sion technique, originally due to Silverstone and MoatsjSq, and introduced to statistical 
mechanics by this author 87], l88|. The general objective of the expansion is to expand all 
functions in terms of spherical coordinates centered at some single point ri. A function of 
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r 23 is naturally expanded in terms of spherical harmonics centered at r 2 or rs; expanding 
a function /(r 23 ) around ri appears unnatural. However, once all functions are expanded 
in terms of spherical harmonics centered at a single point ri together with polynomials in 
the scalar distances ri 2 , ri 3 ,ri„ from particle 1, all angular integrals are integrals over 
products of spherical harmonics, which are basically trivial. By means of the spherical har¬ 
monic technique, a general A^-particle cluster integral can be reduced from 3N to N — 1 


non-trivial integrations 


87l |. While the spherical harmonic expansion can yield an inhnite 


series in the order I of the harmonics, analytic calculations show that the inhnite series typ¬ 
ically converges exponentially rapidly (in 1) to the correct answer 88(]. In equations ll591fTB^ 
matters are simpler. The spherical harmonic expansion here reduces J dr ds to a non-trivial 
two-dimensional integral. 

For the dynamic structure factor of hard spheres, the result 


{g^^\q,0)) ' = 1 + 80-F3O02 + ... 


(166) 


obtains. 

Combining the above equations, for a solution of hard spheres at volume fraction 0 one 
predicts for low q 


Dm = Do{l - 0. 


19.00^ -f- ...) -|- q ■ /S.Dm ■ q. 


(167) 


Some European authors 65|, l68|, l83j] give an alternative form which is approximately 


Dm — Dq{\ -\- 1.560 -|- 0.910^ -|- ...). 


(168) 


Eq. 1 1681 differs from eq. I167l in its 0^ coefficient because the latter equation: neglects all terms 
in V ■ /ijj, assumes that the net solvent volume hux rather than the net total volume hux 
across closed surfaces must vanish, and ignores dynamic friction effects. A demonstration 
comparing our derivation of eq. 11681 with other published results as published in the open 


literature 


89[ | is given in Appendix A. 


H. Dynamic Friction Effects 


Re cog nition of the importance of dynamic friction in diffusion can be traced back to 


Mazo 


90( 1 , who analyzed results of Stigter et ah |9]| on Dg of sodium lauryl sulfate micelles at 
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various surfactant and background electrolyte concentrations. Schurr j92l. l93| gives a related 


treatment of electros 
and the Kirkwood 


;atic effects in Dg. Beginning with the Einstein relation D = ksT/f 


94l | fluctuation-dissipation equation 
1 


/ = 


SksT 


F(0).F(()>, 


(169) 


Mazo showed that / of a micelle has a component due to the micelle-micelle part of F(f); 
his calculations 901] on a simple model for a charged micelle found good agreement with ex¬ 


periment 


9l[ |. these results were anticipated by the corresponding demonstration of Chan¬ 


drasekhar and von Neuma nn [78l| of the existence of dynamic friction in the stellar dynamics 
of star clusters. 

The Kirkwood equation suggests that interparticle forces should also augment the Stokes- 
Law drag coefficient fd, which pertains to straight line unaccelerated motion. In the ap¬ 
proximation Tij ^ 0, this author 95(| found 


sm = 


ksTco 

(27r)3 


' v + rfc)f] -1 


(170) 


gW{k)''“ {k-\ + rk)t 

where 5fd{t) is the time-dependent increment to fd] 'jk = Dmk^ ~ kBT[l — ch{k)]{l — (j))/fo] 
0, c, and fo are the solute volume fraction, concentration, and free-particle drag coefficient, 
respectively, and where h{k) = J dr [g^‘^\r) — 1] exp(ik ■ r). For hard spheres, eq. 11701 shows 
fd = /o(l + ct/0) ct/ = §■ For charged hard spheres with an auxiliary Debye interaction, 
af can become extremely large, as previously shown by Mazo. Similar considerations show 
the contributions of corresponding effects to the mutual diffusion coefficient, the self-diffusion 


coefficient, and the solution viscosity. 


96 


The formalism developed above for Dm, Dg, and Dp leads naturally to the terms 


of eq. 11471 |96| and ADg of eq. 11851 |97j. ADg of eq. 11851 has essentially the same form 
as the Kirkwood equation 11691 differing only in the derivation. In the limit of small q, 
exp(fq ■ Tij) ^ 1 so at small q one has AD^ = ADg. To lowest order in q, dynamic friction 
has the same effect on Dm and Dg. Both AD terms are 0{q^) so neither ADm nor ADg 
vanishes at low q. 

AD may be understood as arising from a caging effect, in which the motions of each 
particle are modihed by interactions with its neighbors. No matter which way a particle 
moves (no matter what the orientation of Ax^i = / \Bi{s)ds), the particle’s neighbors lag 
in responding to the particle’s displacement. On the average, then, for any Ax^j particle 
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i experiences a retarding force that drives the particle back towards its original position. 
However, by diffusion the cage recenters itself on the new location of the diffusing particle, 
the recentering being described expanding the particle’s radial distribution function into its 
sinusoidal and cosinusoidal components, and allowing them to relax via diffusion toward 
their new equilibrium values. 

Physically, dynamic friction acts by dispersing the random Brownian force F Bi on each 
particle over the nearby particles. In the presence of dynamic friction, the Brownian motion 
of each particle is driven by an average over the random forces on a particle and its neigh¬ 
bors. The probability distribution of the averaged force is narrower than the probability 
distribution of the individual F^j, so averaging the random force reduces Dm and Dg. 

The physical origin of dynamic friction is clarihed by considering a particular interaction 
potential between two Brownian particles, namely a covalent bond. A covalent bond serves 
to average the random force on either particle in a dimer over both particles. A random force 
that would have displaced a free particle 1 through Ri will displace (up to corrections arising 
from hydrodynamic interactions) each particle of a covalently bonded pair through Ri/2. 
Superhcially, while the covalent bond retards the motion of either particle, the bond does 
not appear to reduce the mass flux m5x of diffusion; instead of moving a pair of particles 
each of mass m through distances 6xi and 5x2, the random force moves a particle of mass 
2m through a distance |(5xi -|- 5x2). Less superhcially. Brownian forces are random. If the 
free particles would have had random displacements Ri and R 2 , the displacement of either 
particle in a linked pair during the same time interval would be Rd = (Ri -|- R2)/2. The 
random displacements Ri and R 2 are independent random variables; they therefore add 
incoherently. The probability distribution for R^j is thus narrower than the distributions 
for Ri or R 2 , so the random displacements of either particle in a dimer are less than the 
random displacements of a monomer. Forming a dimer reduces the mass hux arising from 
diffusion; in terms of the above treatment, if a monomer of mass m moves a distance 5x the 
covalent bond causes a dimer of mass 2m to move a distance less than 5x/2. 

Several authors have suggested that dynamic friction should affect Dg, without neces¬ 
sarily affecting Dm- For example, Ackerson 9^ applied the Mori formalism to solve the 
Smoluchowski equation. 


9p(r", t) 

m 


N 






(171) 
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obtaining a dynamic friction contribntion to that vanishes if hydrodynamic interactions 
between particles are not included in the calculation. Here p(r^,t) is the iV-particle density 
function, Fj is the direct force on particle j, and is the Smoluchowski two-particle 
diffusion coefficient. In contrast, AD^ is non-zero for hard spheres with no hydrodynamic 
interactions. Implicit to Ackerson’s analysis of the Smoluchowski equation is the assumption 
that the diffusion tensor is a purely hydrodynamic object, all effects of the direct inter¬ 


actions being contained in t 


re sedimentation term Dijf3Fj. However, comparison with the 


Langevin approach shows 6^ that is properly interpreted as a dressed diffusion coeffi¬ 


cient Dfj^ Dij + ADjn, whose value depends in part on direct interactions. The Ackerson 
calculation is correct, except that it does not capture the part of dynamic friction already 
hidden with The Smoluchowski and Langevin equation approaches therefore agree, 

except that the Langevin approach suggests how to calculate ADm but the Smoluchowski 
approach must receive ~ D^j + AD^ from an external source. 

The cage size /, which is determined by the natural range of the interparticle interactions, 
is described by g^‘^\r). For example, two charged spheres of radius a in a solvent of Debye 
length K~^ have / ~ a -|- independent of the sphere concentration. The natural time 
scale on which 6fd attains its long-time value is the time required for the particles to diffuse 
across a cage, which is Td = t^fo/ksT. At times <C Td, Sfd —>■ 0. The longer the range of 
interparticle interactions, the greater will be the time required for 6fd to reach its limiting 
value. Since Z ~ c°, is independent of c. The typical distance between near neighbors 
(that is, the concentration c) determines the strength of the cage (so AD ~ c^), but does 
not directly affect the size of the cage (so Td ~ c°). 

How does one interpret a time-dependent AD^ (and hence a time-dependent hrst cumu- 
lant A'l)? There are three natural time scales involved, namely Td dehned above, the shortest 
time Ti resolved by the correlator, and the time Tq at which the quadratic correction Ct"^ 
to S{q,t) = A — Bt + CD becomes signihcant. Tq is also the time scale on which particle 


displacements become comparable to g so that, e. g., the difference between T, 




Jq)) 

and Tjj(rjj(0)) is signihcant. For meaningful data, one must have ti <C Tq. If S> Tq is 
the longest of these three times (as is found in the high-g limit), the spectrum decays before 
particles encounter dynamic friction; Ki is not modihed by ADm- Conversely, if Ti, 
ADm reaches its plateau value before any data is obtained; Ki contains ADm as a constant. 
If Ti < Td < Tq, the time- dependence of ADm will be visible in the spectrum, leading to a 
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time dependence of the nominal Ki defined above. 

Hess and Klein compute quantities equivalent to AD^ fo r a v ariety of systems, finding 
good agreement with the experiments of Gruener and Lehmann lOOl ] on QELSS of interacting 


polystyrene latex spheres. The model of Hess and Klein[Mi predicts S{q,t) ~ t at long 
time. The experimental S{q,t) ~ for a ~ 1.5 — 1.2, found by this author 30j for 


interacting polystyrene latex spheres in systems of low ionic strength, supports the Hess- 


Klein 


81| model. 


Fluctuation-dissipation calculations related to eq. 11691 predict an increment in solution 
viscosity from particle-particle interactions. In the limit in which eqs. 11471 and 11851 have 
been evaluated, AD^ = ADg, but the viscosity increment is not the same as the increments 
to Dm and Dg. Dynamic friction causes the Stokes-Einstein equation to fail, changes in D 
and r] not being proportional to each other 981. 


I. Microscopic Treatment of Reference Frames 

Reference frame effects are undoubtedly familiar to small children. While using a bathtub, 
most children are all too aware that if they move rapidly from one end of the tub to the 
other, both the water and all the toys floating in the water move rapidly in the opposite 
direction. The objective of this Section is to restate this familiar observation of household 
physics in a slightly more elaborate mathematical form useful for the discussion of light 
scattering methods. The restatement is originally found in Phillies and Wills 8^, which is 
followed closely here. 

The physical basis of reference frame corrections is the fundamental requirement that the 
volume flow of an incompressible solution across a closed boundary must vanish. Container 
walls—through which volume flow is impossible—may be part of the “closed” boundary. For 
a hnite volume of an incompressible solution in a closed container, an equivalent statement of 
the fundamental requirement on volume flow is that volume flow across a plane that bisects 
the volume must vanish. Systems with non-zero volumes of mixing, in which the volume of 
the system is changed by concentration fluctuations, raise fundamental complications not 
included here. 
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The ultimate objective of this discussion is to verify that equation 11631 


dvi [q • T'(ry) ■ qexp(iq ■ Vij) + (l)jH{q)] = 0, 


(172) 


Jv 

is a correct statement of the reference frame correction. In this equation, V is the container 
volume, q denotes the unit vector of q, cj)j is the volume fraction of the single particle 
j which is driving the flow, H{q) is the spatial Fourier transform of the hydrodynamic 
excluded volume of particle j, and T'(rjj) is the exact hydrodynamic interaction tensor 
giving the flow of solution at i arising from motions of particle j. For a solution containing 
N identical solute particles, cf) = N(f)j is the volume fraction of solute in the solution. The 
Oseen tensor T(rjj) is an approximation to T' . The location of the particle driving the flow 
is r^; Vi labels all the locations in V. 

The volume flux has two components. A moving particle transports its own volume across 
any surface it crosses. A moving particle also induces a volume flux in the surrounding 
medium. The volume flux of a particle across a surface S is determined by the particle’s 
velocity and cross-section within S. The volume flux of a particle j across a surface S may 
be written 


Q(5)s-v,. (173) 

Here S' is a plane across the container, s is the normal to S, and Cj{S) is the cross-sectional 
area of particle j contained in the surface S. If particle j is not intersected by S, Cj{S) = 0. 

A moving particle j with velocity Vj and location induces at the point r* in the 
surrounding fluid a fluid flow vp 




) = T'( 


■n. 


■ V 




( 174 ) 


This paper is based on linear hydrodynamics, in which fluid flows due to different particles 
simply add, so to hrst approximation the total flow at rj caused by all the particles in the 
system may be obtained by summing eq. 11741 on j. The flow across S may be obtained by 
integrating Vj ■ s, the solvent flow perpendicular to S over all points of r* in S. 

Choose for a closed surface—across which the volume flux is required to vanish—some of 
the walls of the container, together with a plane S that is perpendicular to the scattering 
vector q. Across this surface, the volume flow is 

■ Vj + [ dS [s • T'ivij) ■ Vj] = 0, (175) 

Js 
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where the f dS includes all points in the plane S. There are no flows across any of the walls 
of the container. 

For q||s, the unit vectors q and s are interchangeable. Within the plane S, exp(iq ■ rj^) 
is a constant. As a result 


/ dviS- T'(rij) • Vj exp(iq • (r^ - rj)) = - dvi C{ri) exp( 2 q • (r^ - ri))s ■ v^, (176) 

Js Js 

Here C'(rj) is 1 or 0 depending on whether or not particle j has part of its volume at the 
point Vi in the surface C. Since particle j is in general of finite extent, the statement that 
particle j has a cross-section that is intersected by S does not imply that the center of mass 
Tj of particle j must lie in S. 

On the Ihs of eq. 11761 (neglecting points near the walls), the only tensors available to 
form T' are I and rr, so by symmetry only the q component of Vj can contribute to the 
integral. Therefore, within the integral on the left hand side one may replace T' ■ Vj with 

T' ■ qq ■ Vj. 

Also, the fluid volume may be decomposed into an (infinite) series of planes parallel to 
S. A sum over the contents of all these planes includes the entire volume of the container. 
Integrating eq. ll76l over all these planes, the rhs of eq. II761 gives the spatial Fourier transform 
(j)jH{q) of the particle volume. Choosing the origin to be at r*, Cj is implicitly a function 
of Tij rather than r*, so it is useful to define 


'V 


dvi exp( 2 q ■ Yij)C{Yi)) = VjH{q), 


(177) 


where H{0) = 1. Defining (f)j = VjjV, and noting that the rhs of eq. 11771 is independent 
of all position coordinates, one may write Vj = f dYi Vj/V] here is a dummy variable of 
integration over a constant. Thus, noting ^aiipianes^^ s = q, one has 


'V 


dr-i q ■ T'(rij) ■ qexp( 2 q ■ rij)(q ■ v^) = -VjH{q){^ ■ v^). 


(178) 


Eq. 11781 is true for arbitrary q ■ Vj, so 


dYi [q ■ T'(rij) ■ + (t)jH{q)\ = 0, 


(179) 


Jv 

completing the desired demonstration. Finally, T ^ T' except at large distances, where 
T' T 0, so eq. 11791 can be subtracted from eq. 11571 to obtain eq. 11641 
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J. Wavevector Dependence of Dm 


Equation 11671 is to be understood as the long-wavelength {q 0) limit of the more 
general forms of eqs. I1551I162I In the earlier equations Ia and Id have only a simple 
dependence, but all other integrals contributing to Dm contain exponential factors exp(—iq- 
r). The exponentials can give Dm a dependence on |g|. If g is small enough, exp(—iq ■ r) 
is approximately 1 or 1 — iq ■ r; the Jj are then all ~ g^, so that Dm ~ g°- However, if g is 
not small (for hard spheres, if ga > 0.1 or so), the exponentials modulate the kernels of the 
integrals; the Jj are no longer simply proportional to g^. Equivalently, at large g the mutual 
diffusion coefficient becomes g-dependent. 

The wavevector dependence of Dm is signihcant when the scattering length q~^ becomes 
comparable with the effective range Og of interparticle forces (this length is, roughly speak¬ 
ing, the distance over which differs appreciably from unity). The third length in 

the problem, namely the mean distance between nearest-neighbor particles (sometimes er¬ 
roneously described as the mean interparticle distance), affects the strength of interparticle 
interactions, but does not directly change the wavevector dependence. Changes in concen¬ 
tration do indirectly change the wavevector dependence by changing g{r). Erroneously? 
A typical distance between two particles is not the mean nearest-neighbor distance, it is 
half the distance across the container. In many systems, gOg -C 1, so Dm is perceptibly 
independent of g. 


Pusey et ah jini| reported the hrst observation of a Dm that is g-dependent due to in¬ 


terparticle forces, namely Dm of charged R17 virus in nearly pure water. At nearly the 


same time, Altenberger and Deutch[l| demonstrated that Dm from the QELSS spectrum 
of a macromolecule solution contains a concentration- and wavevector-dependent correction 
term arising from interparticle interactions. These results have been substantially extended. 
For highly-charged spheres in salt-free water, particles can be sufficiently far apart that hy¬ 
drodynamic interactions are nearly negligible even though electrostatic forces remain large. 
Under these conditions the g-dependence of Dm is predicted to arise primarily from the 
g’^^)(g, 0)“^ term of eq. 11551 so that Dmg^^\q,I^) shoul d be independent of g. This hnal 
prediction is conhrmed by work of Gruener and Lehman 100 1. 


It is sometimes said that a g-dependence of the mutual diffusion coefficient implies that the 
underlying process is not diffusive, by which is meant that Pick’s Law J(r, t) = — ZlmVc(r, t) 
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must be replaced by 


dR Dm{r — R)Vc(R, t). 


J(r,() = 



(180) 


The non-local diffusion tensor -Dm(r “ R-) replaces the local diffusion coefficient Dm- The 
non-local diffusion tensor has a maximum range SR over which it is effective. If one is 
only concerned with diffusion through distances much larger than SR, Dmi^ — R) and Dm 
are indistinguishable. Equivalently, if q~^ S> SR, Dm from QELSS will be perceptibly 
independent of q. 

Eq. [155] becomes especially interesting in the limit q — )■ cxd. In this limit, factors exp(— iq- 
r) and exp(— iq ■ s) oscillate rapidly with respect to the remainder of the integrands of eqs. 
I156II1621 so integrals over these exponentials vanish. For the same reason, at large q the 
distinct terms of g^^\q, 0) tend to zero. The remainder of eq. [T62| vanishes by symmetry, so 


lim Dm = D, 

q—)^oo 


1 + c / * (q ■ bi 2 ■ q) 9 '''qr) + ^ / drds (q ■ bus • q) 9 ' qr, s) +(}■ AD„-q. 

(181) 

Our analysis of ADm reveals that ADm arises from a “caging” or “averaging” effect that 
becomes signihcant when the diffusing particles have diffused through distances comparable 
with the distance over which — 1 is non-zero. In the q ^ oo limit, particle motion is 

only observed at very short times, because at longer times S{q, t) has decayed into the noise 
in the spectrum. At very short times, the particle has only moved through very short dis¬ 
tances. If particles only move over the short distances over which particle motion is observed 
in the q ^ oo limit, ADm is not effective at retarding particle motion; limg_,.oo ADm —t 0. 

Eq. 11801 may be written 


lim Dm = Do{l + fcioo0 + k2oo4>^). (182) 

q—>^oo 


bi 23 is a three-particle term, which first contributes to Dm at the leve 


. The (a/r)' 


approximations to the b give kioo = —1.734 and k 2 oo = 0.873. Batchelor’s 65|] treatment of 
these tensors gets kioo = —1.83. 


K. Self Diffusion Coefficient Dg and Probe Diffusion Coefficient Dp 

Dg and Dp both measure the diffusion of a single particle through a uniform solution. 
They differ in that the self-diffusion coefficient refers to the diffusion of a single particle 
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through a solution of other particles of the same species, while the probe diffusion coefficient 
refers to the diffusion of a given particle through a solution of particles of some other 
species. Dg and Dp are variously measured using FPR, PGSE NMR, PCS, RIGS, QELSS, 
and macroscopic tracer diffusion techniques, with some techniques only measuring Dg and 
others only measuring Dp. For these experiments, and taking af = 1 and thus gig{q, 0) = N 
for mathematical clarity, 

r) / 1 ^ 

Ds = * Vj(s)Vi((+ t) + ...) 

(183) 

Replacing Vj(t) with its Brownian and direct components, and applying eqs. 1130111441 

TV 1 ^ \ 

I + ^ 1 ■ q + q ■ ^Dg ■ q (184) 




where 

(185) 

As q ^ oo, ADg vanishes, because particles have not yet moved a distance comparable to 
the range of g^‘^\r). Eqs. 11841 and 11811 then reveal that in the high-g limit the concentration 
dependences of Dg and Dm are the same. Furthermore, to C>(g^), ADm includes only the 
self {i = j) terms of AH’", which are the same as the individual terms of ADg (eq. I185p . so 
to 0{q‘^) the dynamic friction contributions to Dg and Dm are equal. 

Optical probe experiments are based on ternary solvent: matrix: probe systems. In 
applying QELSS to such a system, in many cases one works in the limit that the probe 
species, while dilute, completely dominates scattering, while the matrix species scatters 
negligible amounts of light, even if it is concentrated. In other cases, measurement of the 
QELSS spectrum of the probe:matrix solution and of the probe-free matrix solution, followed 
by subtraction of the latter from the former at the field correlation function level, removes 
the matrix scattering from the scattering by the mixture, permitting successful isolation of 
the probe spectrum |l3l|. 

The above formalism can evaluate Dp if some minor adaptations are made. Namely, in 
eq. O for a probe experiment the scattering lengths cxj assumes two values: CTj = 1 for probe 
molecules, and Uj ~ 0 for matrix molecules. Sums over particles may include either matrix 
or probe molecules. The Nm matrix molecules are labelled {1, 2,..., Nm}, while the Np probe 


N 


ADg = ^ \ X] / ds (vBi(s)vDj(t) + V£,i(s)vBi(f)) 


, i=l 
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molecules are labelled Nm + 1, Nm + 2, Nm + Np^ the total number of molecules being 
Nt = Nm + Np. The matrix and probe concentrations are cm = Nm/V and Cp = Np/V, 
respectively; for dilute-probe results only terms of lowest order in Np or Cp are retained. Since 
the matrix and probe species are not necessarily the same, the hydrodynamic interaction 
tensors b and T and the radial distribution functions may have different values for 
probe-probe, probe-matrix, and matrix-matrix molecular pairs. 

Under these conditions, eq. 11511 gives us the probe diffusion coefficient 

/ Np Nt Np 

aiaj exp(iq ■ rij)iqVi : [Dji] + ^ crjCTj-q ■ Aj ■ qexp[iq ■ ri^]) 

^ ^ 1 \ ij=i ij=i 


Np 

+ CTiCTj-q ■ AAj ■ q 

Terms with i,j > Np vanish because cTj, aj ~ 0 for matrix molecules. While matrix molecules 
are optically inert, they are still hydrodynamically active, so neither nor Tj; vanishes if 
I refers to a matrix molecule. 

The signihcance of eq. 11861 becomes clearer if YA- ■ is resolved into self and distinct com¬ 
ponents, and limited to terms referring to no more than two particles, viz: 

Np / Nrp 

Y ■ I ■ q - 5^(q ■ biz ■ q) + zqVi : [I + b^] 

i=l V 1=1 m=l 


D = 



N’j' 


AD.. 


+iq^ Vz : [Tiz]) - q - ■ q 


1=1 


+ 


Np 

(TiCTj 

i^j=l 

In this equation, is the single-particle diffusion coefficient. Terms in cx^iqV vanish by 
symmetry. Terms in exp{iq ■ r^j) are only non-zero if two probe particles are close 
enough to interact (close enough that their ^ 1). If the probes are dilute, this event 
rarely happens, so the aiaj terms are negligible with respect to the af terms. With these 
reductions. Dp becomes 



(exp(iq ■ rij)(-q ' T^ji ' q + ^qV* : [T^z] -h iqV^ : [I])) 




. (187) 


D - 

" «V‘>(9.0) 


■ Np 

E 

2=1 


(T, 


N'T 


AD.. 


-q • I ■ q - ^ q ■ (b^z) ■ q - q ■ ■ q 


1=1 


(188) 


Through the concentration effects on Dp, Dg, and Dmi^q —>■ oo) are now seen to be 

identical. 
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V. OTHER APPROACHES 


A. Coupling of Concentration and Energy-Density Fluctnations 

A pure simple liquid such as water has an intrinsic spectrum arising from scatter¬ 
ing of light by propagating pressure fluctuations and non-propagating energy-density 
fluctuations Scattering from sound waves creates Brillouin peaks, whose centers are 
displaced from the incident light frequency. Scatteri ng fro m energy fluctuations (sometimes 
described as “entropy” or “temperature” fluctuations [l02j) leads to a central “Rayleigh” line 
with width proportional to the thermal conductivity S. If a solute is added, the spectrum of 
the pure fluid gains a further line arising from concentration fluctuations; this mass-diffusive 
line is the spectrum treated above. It is almost always true that scattering by the pure liq¬ 
uid is much weaker than scattering by concentration fluctuations, so that concentration 
scattering dominates the spectrum. 

In general, concentration and energy fluctuations are coupled, as by the Soret effect (a 
temperature gradient driving a mass flux) and the Dufour effect (a concentration gradient 
driving heat flow). Correlations in pressure, energy-density, and concentra tion fluctuations 


are treated by Mountain and Deutch, and by Phillies and Kivelson 


103, 


104| . Cross-coupling 


of energy and concentration diffusion modihes the relaxation rate of heat- and mass- dif¬ 
fusive modes, and creates an additional spectral line (of integrated intensity zero) due to 
static cross-correlations between energy and concentration fluctuations. If 
and if solvent scattering is weak (conditions readily satished by macromolecule solutions), 
the spectrum decouples, so that the relaxation of concentration fluctuations is determined 
entirely by Dm- 

B. Smoluchowski and Mori-Zwanzig Formalisms 


The purpose of this short section is to put in one place discussions on the Smoluchowski 
diffusion-sedimentation equation and outcomes from the Mori-Zwanzig formalism. The 
Smoluchowski equation 


9p(r", t) 

m 


N 






'dvi 


(189) 
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connects the time dependence of the iV-particle probability distribution function 
to a set of external forces Fj on the particles and to a diffusion coefficient The 

diffusion coefficient is usually taken to be the hydrodynamic diffusion coefficient Dij = 
ksTpij, The diffusion-sediment at ion equation was obtained for systems in which the forces 
are imposed externally, so that there are no correlations between the forces and the Brownian 
displacements. In this case the dynamic friction term vanishes, an outcome that is known 


experimentally 90|, l9l| to be incorrect for diffusing systems. In order apply the Smoluchowski 
equation to a diffusing system, for starters one needs to reinterpret to included dynamic 
friction terms. Such an approach is not considered further here. 

In the Mori-Zwanzig formalism, transport coefficients are taken to arise from force-force 
correlation functions (memory kernels) that describe the temporal evolution of the projected 
force. In conventional Mori-Zwanzig calculations, one sets up the transport equations, but at 
the key point the memory kernels for the Mori-Zwanzig projected forces are simply identihed 
as the appropriate transport coefficients. The Mori-Zwanzig formalism yields a transport 
coefficient, but the meaning of that coefficient is dehned by the Mori-Zwanzig equation. The 
transport coefficient’s value is obtained by htting experimental measurements to solutions 
to the Mori-Zwanzig solutions. The symbols used for Mor-Zwanzig transport coefficients are 
the same as the symbols used in classical measurements, but that similarity may be shallow. 

One can simply say that the Mori memory kernel corresponding to mutual diffusion is 
for some Vij, but having done so one cannot also claim that Vij is ksTfiij with pij being the 
hydrodynamic p. That would be overdehning Any expression other than ksTpij that 
has the correct dimensions and g-dependence is also consistent with the unevaluated Mori 
Kernel. In particular, the rhs of eq. 11461 including the dynamic friction terms of eq. 11481 and 
the terms in is the correct evaluation for the lead terms in q of the Mori kernel. 

Alternatives to simple identihcation of particular memory kernels as particular trans¬ 
port coefficients do exi st. A Mori kernel has been extracted from computer simulations 
on a simple system 105(]. The memory kernel for the projected forces, including the use of 
the proj ected time evolution operator, has been evaluated analytically by Nasto in a spe¬ 


cial case 


1061]. Direct evaluations of the Mori memory kernel to determine hydrodynamic 


transport coefficients corresponding to Dfj^ have not been made. 

Finally, we note the notion of using the Mori formalism to generate solutions of the 
Smoluchowski equation. This approach is not justihable in terms of derivations of the 
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Mori formalism. The Smoluchowski equation has a diffusion term, driving concentration 
fluctuations to decay toward zero, but unlike the Langevin equation has no random force 
term . As a results, solutions of the Smoluchowski equation are not stationary in time; they 
show, e.g., ag{t) decaying to zero. In contrast, derivations of the Mori equation explicitly or 
implicitly (for example, by assuming that certain Laplace transforms exist) assume that the 
Mori equation solutions are stationary in time. In consequence, one can formally insert the 
Smoluchowski equation into the Mori formalism, but the meaning of the outcome of such a 
process is unclear. 


VI. DISCUSSION 


A. Implications for QELSS Measurements 


The diffusion coefficients measured by QELSS, FPR, PCS, RIGS, PGSE NMR, inelastic 
neutron scattering, and other techniques all depend on the solute concentration. The treat¬ 
ment above shows how these dependences may be calculated. The Stokes-Einstein equation 

keT 


D = 


6n'r]a 


(190) 


is sometimes used to interpret Dm from QELSS measurements. For dilute neutral spheres, 
radius a, diffusing through a simple solvent having a modest viscosity rj not too much larger 
than the viscosity of water, this equation can be appropriate. The equation is assuredly 
invalid for macromolecules diffusing at elevated concentrations. Equivalently, Dm, ?7, T, 
and eq. 11901 cannot be combined to calculate the hydrodynamic radius of non-dilute dif¬ 
fusing spheres^ _ Dm in non-dilute solutions is nearly certainly not independent of solute 

concentration lOSj, though over limited ranges of c the dependence may be small. 


We considered in detail the mutual and self diffusion coefficients of neutral hard spheres. 
The lead term in the dependence of Dm on solute volume fraction (j) (eq. I167p is small. There 
are large hydrodynamic and thermodynamic contributions to dDm/dc, but for neutral hard 
spheres at concentrations that are not too large these contributions almost cancel. Changing 
the intermacromolecular interactions, for example by charging the spheres, will change the 
hydrodynamic and thermodynamic contributions to dDm/dc. There is no reason to expect 
the changes in the different contributions to cancel. For systems that are not neutral hard 
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spheres Dm may depend appreciable on concentration. Indeed, the dependence founc 


charged macromolecules can be much larger than the dependence for neutral spheres 
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for 


79|. 


Particle sizes in concentrated solution have sometimes been estimated from QELSS data. 
The particle shapes and interactions must be known. The inferred radii depend on the 
detailed theoretical model in use. 


59 


60| . Measurements of probe diffusion by polystyrene 


latex spheres of different sizes diffusing through micelles, followed by interpretation applying 
the above hydrodynamic models, have allowed determination of the (substantial) water 
content of Triton X-100 and other micelles, jl^ 


van Megen and Underwood 


109[ | report a light scattering method that can determine 


the distinct part of the held correlation function Their approach is based on 

continuous variation in the solvent index of refraction, coupled with the use of two colloidal 
species having the same size and surface properties, but different cores and hence indices 
of refraction. The same physical approach is found in inelastic neutron scattering, using 
deuterated and hydrogenated macromolecular species and a series of mixed deuterated and 
hydrogenated solvents. 


B. Comparison of Dm and Dg with Experiment 

This section treats experimental data suita ble fo r testing the aforementioned theoretical 
models. There is also an extensive literature lin| on diffusion, crystallization, and glass 
formation in concentrated sphere suspensions. However, the concentration regime in which 
spheres vitrify, and the concentration regime in which the above calculations are likely to 
be valid without extension to higher order in 0, are not overlapping. 

die best study of Dm at low g in a hard-sphere system appears to be that of Mos et 


al. 


lll| . who used homo dyne coincidence spectroscopy (HCS) to measure Dm- Homodyne 


coincidence spectroscopy 


112| is a two-detector QELSS experiment, in which a sample is 


simultaneously illuminated with two incident laser beams, the scattered light is collected by 
two detectors placed on opposite sides of the scattering volume, and the inten sity- intensity 
time cross-correlation function is measured. As was hrst shown by this author 113| , the ho¬ 
modyne coincidence spectrometer differs from a conventional one-beam one-detector QELSS 
instrument in that an HCS spectrometer is substantially immune to multiple scattering ar¬ 
tifacts. While multiply-scattered light reaches each HCS detector, only the single-scattered 
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light has cross-correlations from one detector to the other, so the only time-dependent com¬ 
ponent of an HCS spectrnm is that dne to single scattering. 


Mos et al. 
stearylatedpjj 


ll| stndied colloidal silica spheres, prepared by the Stoeber process |ll4j| . 


, and snspended with sonication in xylene and tolnene at concentrations of 


0 - 180 gjL. Over this density range, exhibits a linear decline with increasing c, i.e.. 


Dm = Do{l a0). 


(191) 


Using the estimated density of 1.75g/cm^, Mos et al.’s data for xylene solntions shows 
a = —0.86; in tolnene solntions, a = —1.2 was fonnd. Here Do is the zero-concentration 
limit of Dm and (j) is the sphere concentration in volnme fraction nnits. The total change 
in Dm over the observed concentration range is ca. 10% of Do- Inferring from the data on 
Dm and & 3-5% nncertainty in Mos’s individnal measnrements of Dm, the nncertainty in 
a mnst be ronghly ±0.1. The theoretical a of eq. 11671 is -0.9, in excellent agreement with 
these experiments. 

The largest nncertainty in the experiment is the identihcation of stearylated silica par¬ 
ticles, th ~ 370A, with a Cig chain coating, as hard spheres. If the particles had weak 
attractive or repnlsive interactions, in addition to their hard-sphere interaction, the con¬ 
centration dependence of Dm would be expected to be different. An attractive interaction 


usually makes a more negative, while a repulsive interaction usual 
tive. Thus, if one wished to believe the Batchelor-Felderhof solution 


y makes a more posi- 


65 


^1 of the diffusion 


problem (which leads to a = ±1.56 for hard spheres), one could claim that stearylated silica 
particles attract each other weakly, thereby reducing their a from the Batchelor-Felderhof 
value to the value found experimentally. Mos, et al., in fact make this i nterp retation for 


their spheres. Caution is needed to avoid circular arguments. Mos, et al.. 


lll| began with 


the assumption that the Batchelor-Felderhof calculation is correct, and therefore sensibly 
inferred from their observed negative a that their spheres must attract each other weakly. 
However, they do not in their paper adduce other evidence that the spheres had an attractive 
potential, in addition to the hard sphere potential. If eq. 11671 were correct, the a found for 
stearylated silica spheres would have the value appropriate for hard spheres, consistent with 
the physical expectation that these uncharged spherical objects that form stable solutions 
in non-polar solvents shoul d be close to true hard spheres. 


Kops-Werkhoven et al. 


116| used QELSS on stearylated silica spheres in cyclohexane 
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to determine Dm at low g, finding a = +1.56, in agreement with the Batchelor-Feld erho 


theory 


65 


80| . The spheres used by Kops-Werkhoven et ah jll6| and by Mos et ah 


111 


are chemically the same, so in hrst approximation both groups either did or did not study 
hard spheres. The sphere-sphere interactions could differ slightly, because the two studies 
used different organic solvents. Kops-Werkhoven et ah also employed a contrast-matching 
technique to estimate as in 

Dg = DoiX T + ...), (192) 

reporting as = —2.7 ± 0.5, which is not in good agreement with the theoretical (hydro- 
dynamic) estimate as ~ —1.83. The agreement of the experimental ag might improve if 
dynamic friction at some level were included in the theoretical estimate for a g. 

For the stearylated silica spheres studied by both groups, Mos et al. jlll| demonstrated 
in a non-polar solvent having an only slightl y les s favorable index-of-refraction match than 


the solvent used by Kops-Werkhoven et al. 


116l | {6n ^ 0.05 in toluene, vs. 6n ^ 0.015 


in cyclohexane, based on cycloheptane being an index-matching fluid for these spheres' 


that artifactual positive values of am can arise from multiple scattering. Mos et al. 


ml 


further showed that multiple scattering is important in spheredoluene mixtures. Mos et ah’s 
measurements on silicamonpolar solvent mixtures, which employed HCS (a techn ique that 
is immune to multiple scattering artifacts), are therefore preferable to ref. 116| ’s QELSS 
measurements. The difference in am between the two sets of measurements is small but 
signihcant, implying weak multiple scattering. 

Qiu, et al.. |ll7l | report Dm of polystyrene latices in water at elevated concentrations (up to 
0 0.45), using diffusing wave spectroscopy (DWS) to study highly turbid solutions. DWS 

is a multiple-scattering technique that measures Dm at very large effective q. However, the 
underlying theory for DWS assumes that the displacement probability distribution P{x, r) 
for the moving particles is a Gaussian, which need not be the casej^, so some caution is 
needed in interpreting these results. The lattices in the study had th of 0.206 fi and 0.456/i; 
DWS determined Dm for particle motions < 50A. The latex particles are here diffusing 
dista nces <C r^, so in this experiment t and the dynamic friction term ADg ~ 0. Qiu 


et al. 


117 | report Dm = Do(l —(1.86+0.07)0), i.e. ^ —1.86+0.07, in good agreement with 


Batchelor’s 65| theoretical estimate fcioo ~ — 1.83 and with our estimate fcioo ~ —1.734. In 
contrast, i n Ko ps-Werkhofen’s determinations jll6 | of particles diffused distances ~ th, 


so in refs. 


116| ADg could have been substantial, perhaps explaining why Kops-Werkhofen, 
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et al., found a more negative value for ag than did Qiu, et al. 

The micr osco pic and continuum models have been extended by Borsali and 
collaborators 118| to systems containing a solvent and two physically distinct random-coil 


polymer species. Benmouna, Borsali, and collaborators 


118| compute spectra of ternary sys¬ 


tems in which one or both components is concentrated and in which one or both components 
scatter light. Their polymer models are somewhat remote from the emphases of this review, 
but the agreement between their predictions and their experiments is excellent, supporting 
the belief that the above results are fundamentally correct. 


Appendix A: Other Methods for Calculating Dm 

The purpose of this Appendix is to compare possible methods 
vant methods include those proposed by Batchelor Felderhof 


or computing D„ 
], and this author 


hele- 


63 


82|. 


These calculations are all correct, but do not all correspond to the same experiment or 


physical quantity. While I have previously shown 8^ why refs. 


65l and 


83 


do not hnd the 


concentration dependence of Dm as measured by light scattering, this demonstration has 
not been widely understood. For example there has sometimes been a misapprehension that 


refs. 


65 


m 


82 


and 


63l only disagree about dynamic friction effects. In fact these works 


disagree as to the correct form for the concentration dependence of Dm, even if dynamic 
friction is neglected. It will here be shown that each reference gives an algebraically correct 
calculation of some diffusion coefficient, but that only refs, 
diffusion coefficient Dm that is measured by QELSS. 


63 


and 


82 


calculate the mutual 


Experimentally, light scattering spectroscopy is directly sensitive to collective coordinates 
determined by particle positions. Specihcally, the instantaneous scattered held is propor¬ 
tional to the instantaneous value of the spatial Fourier component ek{t) of the local index 
of refraction, where ek{t) is in turn proportional to the spatial Fourier component ak{t) 
of the concentration of scattering particles. Here k is the scattering vector selected by the 
source and detector positions. In light scattering spectroscopy. Dm is obtained from the 
temporal evolution of akit), namely 

limt^o (^) (a_fc(0)afc(t)) 


Dm 


In a real experiment, eq. 


-P(afc(0)afc(0)) 

is typically applied via a cumulant expansion of S{k,t). 


(Al) 
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In contrast to eq. lATl Batchelor}^ obtains Dm by evaluating the flux of particles due to 
an applied steady (thermodynamic) force, the flux being related to the diffusion coefficient 
by 

J = —DmVc. (A2) 


In ref. 65| Dm was obtained by generalizing the Einstein expression 


D = kBT{h) 


(A3) 


for Dm, b being the mobility tensor. In Einstein’s original analysis, b was a constant. 
In Batchelor’s generalization of Einstein’s analysis, b is given by a microscopic expression 
which depends on the relative positions of the particles in the system. Dm is a macroscopic 
quantity which does not depend explicitly on particle positions. To obtain a macroscopic 
Dm from a microscopic b, Batchelor took an ensemble average of b over possible particle 
configurations. Batchelor’s generalization of the Einstein form for Dm was therefore 


Dm = kBT{h). 


(A4) 


While the calculation of Dm as kBT{h) is mathematically correct, QELSS does not deter¬ 
mine Dm by measuring a flux and a concentration gradient, and taking a ratio. Instead, Dm 
from a QELSS measurement is obtained from the first cumulant in an expansion of t), 

i. e., from the time rate of change of a concentration. While eq. IA4I might happen to give 
the value of Dm that is obtained by light scattering, eqs. lAll and IA2I are not the same. Any 
disagreement between results obtained from these equations is properly resolved in favor of 
the equation that correctly models the quantity that is measured experimentally by QELSS, 
this being eq. lAl 


References (82j and 83(] both made assumptions equivalent to assuming that particle 
motions are correctly described by the Smoluchowski sedimentation equation 

(ic(r, t) 


dt 


= Sc(r,t), 


(AS) 


where the Smoluchowski operator is 


S = V-(D-V + D-F/A;bT). (A6) 

Here D is the diffusion tensor (as distinct from the inferred diffusion coefficient Dm) and F is 
the applied force on a particle. The dynamic friction effects discussed above, which are not 
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at issue here, may be incorporated in D, so that eq. lASl is formally valid no matter whether 
or not dynamic friction effects are present. Eq uat ion I A6 1 manifest Iv includes terms in V-[D], 
which are non-vanishing if D depends on position, or if D depends on concentration and c 
depends on position. 

Just as eq. IA3I was applied to the problem by interpret!^ b as a microscopic mobility 
tensor, so also eqs. lASI and lAhl were applied in refs. |82| and |83| by giving D, F, and c(r, t) a 
microscopic interpretation. In particular, c{r,t) was replaced with ak{t) and its microscopic 
representation in terms of particle positions, namely 


N 

= X]exp(ik ■Tj), 


(A7) 


the sum being over all N particles in the system. is then obtained from some ensemble 
average over the short-time limit of S. 

How do the various applications of eq. IA6I differ from each other? Felderhof’s calculation 


80l | of Dm sets the Smoluchowski equation in the form 
da{ri,t) 


dt 


= U0&/V1 


Via{vi,t) + c(3 I ViH 5 f^^^a(r 2 ,t)dr 2 


+ cWi ■ J hg^‘^^dr 2 ■ Via(ri, t) CqVi ' J t)dr 2 + O(c^). (A8) 

where a(ri,f) is the local density at ri, V is the interparticle potential, = g^‘^\r 2 — ri) 
is the equilibrium pair distribution function, and V* is the gradient with respect to particle 
i, the particles being the particle of interest 1 and a neighboring particle 2. As interacting 
particles must be relatively close together, the concentration gradients at their locations 
should be roughly equal, i. e. 

Via(ri,t) = V2a(r2,t) (A9) 

for a pair of interacting particles. (This is a small-k approximation.) Rearrangement of eq. 
gives 

da{ri, t) 


dt 


as evaluated in ref. 


= I^oVi 


Via(ri,t)-t-c/3 / dr2ViV{Ti2)g^‘^\ri2)a{r2,t) 


+ c (h + T) : Vla{ri,t)g^‘^\Ti 2 )dr 2 


83| . plus terms such as 

c [ Vi- [T] 5 ((^)(ri 2 ) ■ V 2 a(r 2 ,t)dr 2 


(AlO) 


(All) 












which vanish because f V-Trfr is odd in r. Reference 83|] thus finds that terms in V- (T + b) 
do not contribute to Dm- 

References [8^, |82|, and |63| (that is, equations lAlOl and IlSUp disagree because they 


obtain Dm by taking ensemble averages over different functions. In eq. 11501 D was obtained 
from an average having the general form 


(a_fc(0)S'afc(t)) (A12) 

in which is a time evolution operator related to S, while the are statistical weights, 
appearing because the contribution of a particular particle configuration to the measured 
Dm is weighted by the light scattering intensity contributed by that particle configuration. 

In contrast to eq. IA121 eqs. lAlOl and lAllI obtain Dm from an ensemble average over the 
algebraic kernel of S, with no factors of included within the average. The kernel of the 
microscopic S does give the concentration changes to be expected from a given microscopic 
particle configuration, so (S) does give an average rate of change for akit); also, (S) has no 
divergence (V ■ T) terms. However, QELSS obtains the light-scattering-intensity-weighted 
average (the z-average) temporal evolution of ak(t), not the unweighted average. States 
which scatter no light make no contribution to the observed temporal evolution of ak{t). 
The z-weighting in eq. IA12I arises from factors a_fc(0) and ak{t). To obtain a properly z- 
weighted average, factors of ak must be included in the ensemble average. Including these 
factors replaces, e. g., eq. lAlll bv 

c j dr 2 Vi-{D^-^^^)V-[T], (A13) 

a term previously seen in eq. 11501 Unlike the term eq. lAlll the term kV ■ Texp(ik ■ r) of 
eq. IA13I includes parts that are not odd in r and do not vanish on performing the integral. 
QELSS measures {a{0)da{t)/dt), not J, so the appropriate microscopic average for Dm is 
that of eq. 11501 not the form of eqs. IA3I or lAlOl but IA3I or lAlOl do represent diffusion 
coefficients which some experiment measures. 


Appendix B: A Partial Bibliography—Theory of Particle Diffusion 

The following is an incomplete bibliography of pre-1990 theoretical papers on the diffusion 
of mesoscopic particles (colloids, proteins, micelles) at elevated concentrations. I have added 
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a few more recent papers. I make no claim of completeness, though I have tried to be sure 
that most major lines of work are represented. Additions and corrections are welcome. 
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